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ABSTRACT 

We formulate a general method for perturbative evaluations of statistics of smoothed cosmic fields, and 
provide useful formulas in application of the perturbation theory to various statistics. This formalism 
is an extensive generalization of the method used by Matsubara (1994) who derived a weakly nonlinear 
formula of the genus statistic in a 3D density field. After describing the general method, we apply 
the formalism to a series of statistics, including genus statistics, level-crossing statistics, Minkowski 
functionals, and a density extrema statistic, regardless of the dimensions in which each statistic is defined. 
q \ The relation between the Minkowski functionals and other geometrical statistics is clarified. These 

statistics can be applied to several cosmic fields, including 3D density field, 3D velocity field, 2D projected 
density field, and so forth. The results are detailed for second order theory of the formalism. The effect 
t"**- , of the bias is discussed. The statistics of smoothed cosmic fields as functions of rescaled threshold 

by volume-fraction are discussed in the framework of second-order perturbation theory. In CDM-like 
models, their functional deviations from linear predictions plotted against the rescaled threshold are 
generally much smaller than that plotted against the direct threshold. There is still slight meat-ball shift 
against rescaled threshold, which is characterized by asymmetry in depths of troughs in the genus curve. 
\ A theory-motivated asymmetry factor in genus curve is proposed. 

Subject headings: cosmology: theory — large-scale structure of universe — methods: statistical 

O ' 
O 
O 

Any conceivable theory of structure formation in the universe predicts statistical properties of observable quantities. 
O ■ Therefore, the analysis of the present inhomogeneity of the universe is inevitably statistical. However, what kind of 
I ' statistics is useful is not obvious, since we can invent infinite number of statistics to be analyzed. Whether we can adopt 
better statistical descriptions is of great importance in this sense. Each statistic has both advantages and disadvantages. 
The power spectrum, for example, can fully characterize the random Gaussian fields, while it does not contain any 
^ ' information on the non-Gaussianity, which contains significant information in the gravitational instability theory. Among 
various statistical quantities, there is a promising class of statistics which utilizes smoothed cosmic fields. The smoothed 
field has less noisy property than the actual (unsmoothed, or raw) cosmic fields, such as galaxy distributions, temperature 
fluctuations of cosmic microwave background (CMB), shear fields of the gravitational lensing, and so forth. Perhaps the 
; i ■ simplest example of such statistics is the variance (5r ) of smoothed density contrast <5r, which is a function of smoothing 
f3 i length R. Similarly, higher-order cumulants (5r N ) c (N — 3,4,...) are also simple statistics. The density probability 
distribution function (PDF) P(5r), which in principle can be constructed from the hierarchy of cumulants (e.g., Balian 
& Schaeffer 1989), is an another example of popular statistics of smoothed cosmic fields. 

Rather recently, more complex statistics of smoothed cosmic fields have become popular in cosmology, such as the 
genus statistic (Gott, Melott & Dickinson 1986), density peak statistics (Bardeen et al. 1986), area, length, level-crossing 
statistics (Ryden 1988a), Minkowski functionals (Schmalzing & Buchert 1997), etc. These statistics provide assuring 
characterizations of the clustering pattern that can not be perceived only by the hierarchy of cumulants or by the PDF. 
The genus statistic is a powerful measure of the morphology in the 3-dimensional (Gott et al. 1989; Moore et al. 1992; Park, 
Gott & da Costa 1992; Beaky, Scherrer & Villumsen 1992; Rhoads, Gott & Postman 1994; Vogeley et al. 1994; Protogeros 
& Weinberg 1997; Sahni, Sathyaprakash & Shandarin 1997; Canavezes et al. 1998; Springel et al. 1998; Colley et al. 
2000), and 2-dimensional (Melott et al. 1989; Gott, Mao, Park & Lahav 1992; Coles & Plionis 1991; Plionis, Valdarnini 
& Coles 1992; Davies & Coles 1993; Coles et al. 1993; Colley 1997) clustering of galaxies and clusters, and also of the 
pattern of temperature fluctuations of CMB radiation (Bond & Efstathiou 1987; Torres 1994; Smoot et al. 1994; Torres, 
Cayon, Martinez- Gonzalez & Sanz 1995; Park et al. 1998). The peak statistics of the 3-dimensional density field are 
frequently used in connection with the statistics of the collapsing object (Kaiser 1984; Mann, Heavens & Peacock 1993; 
Croft & Efstathiou 1994; Watanabe, Matsubara & Suto 1994; Cen 1998; Gabrielli, Labini & Durrer 2000), while the peak 
statistics of the CMB fluctuations (Sazhin 1985; Bond & Efstathiou 1987; Coles & Barrow 1987; Vittorio & Juszkiewicz 
1987; Kogut et al. 1995; Cayon & Smoot 1995; Fabbri & Torres 1996; Heavens & Sheth 1999) and of the weak lensing 
fields (van Waerbeke 2000; Jain & Van Waerbeke 2000) are suitable for constraining cosmological models. The area, 
length, and level-crossing statistics directly quantify the amount of contour surfaces (Ryden 1988b; Ryden et al. 1989; 
Torres 1994). The Minkowski functionals (Mecke, Buchert & Wagner 1994; Sahni, Sathyaprakash, & Shandarin 1998; 
Sathyaprakash, Sahni, & Shandarin 1998; Kerscher, Schmalzing, Buchert & Wagner 1998; Bharadwaj et al. 2000), which 
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are closely related to the above statistics, are also applied to smoothed cosmic fields (Winitzki & Kosowsky 1997; Naselsky 
& Novikov 1998; Schmalzing & Gorski 1998; Schmalzing et al. 1999; Schmalzing & Diaferio 2000). These statistics of 
smoothed density fields are considered as powerful descriptors of the statistical information of the universe. It is therefore 
essential to establish theoretical predictions of the behavior of such statistics so that we can effectively and ideally analyze 
the data of our universe. 

Recent developments of the perturbation theory (for review, see Bernardeau, Colombi, Gaztanaga, & Scoccimarro 2002) 
in calculating the variance, the cumulants, and the PDF are remarkable. The perturbation theory becomes more and more 
useful because the recent developments of observations enable us to have widely covered sample volume of the universe, 
which can minimize undesirable strongly nonlinear effects which we do not understand well. The direct comparison of the 
theoretical predictions of the perturbation theory with the actual data is a promising field of research in this sense. In the 
case of top-hat smoothing function, Juszkicwicz, Bouchet & Colombi (1993) and Bernardeau (1994a) used the tree-level 
(i.e., lowest-order) perturbation theory to obtain the third- and forth-order cumulants, i.e., the skewness and the kurtosis. 
The same quantities with Gaussian smoothing are calculated by Goroff et al. (1986), Matsubara (1994) and Lokas et al. 
(1995). Quite cleverly, Bernardeau (1994b) took advantage of special properties of the top-hat smoothing function, and 
succeeded to obtain full hierarchy of higher order moments in the tree-level perturbation theory. He also obtained the 
PDF from this hierarchy of cumulants, which remarkably describes the nonlinear behavior of the gravitational instability 
in numerical simulations. Beyond the tree-level calculation, the perturbation theory with loop-corrections has been also 
developed (Juszkiewicz 1981; Vishniac 1983; Juszkiewicz, Sonoda & Barrow 1984; Coles 1990; Suto & Sasaki 1991; Makino, 
Sasaki & Suto 1992; Jain & Bertschingcr 1994; Baugh & Efstathiou 1994; Scoccimarro & Frieman 1996a, b). Because of 
the simplicity of the statistics, calculating the variance, the cumulants and the PDF were primarily the playground of 
perturbation-theorists. The evaluation of other statistics by the perturbation theory is less trivial. A quite useful technique 
is the Edgeworth expansion, which was first applied to cosmology in calculating the PDF from seeds perturbations by 
Schcrrer & Bertschinger (1991), and from non-linear perturbation theory by Juszkiewicz et al. (1995), and Bernardeau & 
Kofman (1995). The analytic expression of the genus statistics in the perturbation theory is derived by Matsubara (1994), 
whose technique corresponds to multivariate version of the Edgeworth expansion. In some literatures the Edgeworth 
expansion is used in connecting the statistics and dynamics of the universe (Chodorowski & Lokas 1997; Lokas 1998; 
Taruya & Soda 1999). The purpose of this paper is to give a comprehensive description of the formalism by which the 
perturbative evaluation is possible for wide range of non-trivial statistics of smoothed cosmic fields in general. 

Since the number of spatial dimensions of our universe is three, the statistics of large-scale cosmic fields are defined in 
cither one-, two-, or three-dimensional space. For example, the space of the density field in a redshift map of galaxies or 
quasars is three-dimensional. A projected galaxy map on the sky, a shear field of gravitational lensing, and temperature 
fluctuations of CMB on the sky are fields in two-dimensional space. The absorption lines of quasar spectra and pencil- 
beam surveys of galaxies define fields in one-dimensional space. Thus it is useful to develop our statistical method in 
d-dimensions for generality. As illustrative examples of applications of our method, we calculate level-crossing statistic 
(or equivalently length and area statistics), genus statistics, density-extrema statistics and Minkowski functional. As 
cosmic fields, we consider density and velocity fields in three dimensions, the projected density field of galaxies in two 
dimensions. The basic formalism and results of the second order theory are presented in this paper. We will give results 
of the third order theory in a future paper, which are technically more involved. 

The primary purpose of the present paper is to fully describe the basic formalism for the statistics of smoothed 
cosmic fields in perturbation theory. Applications of the second order theory to popular statistics and cosmic fields are 
systematically presented. Thus, this paper is in some way a mixture of the new results and comprehensive review of the 
old results. The major new results in this paper are: 

• Explicit formulas of the lowest non-Gaussian correction in various smoothed cosmic fields in <i-dimensions, 
including all the Minkowski functionals. They are presented as functions of both density threshold and rescaled 
threshold by volume-fraction. 

• Relations among several statistics of the smoothed field are clarified. 

• Derivatives of skewness parameter of the velocity field and projected 2D field are derived in second order 
perturbation theory of gravitational instability theory. Especially, skewness parameters with Gaussian smoothing 
are detailed. These quantities are particularly important to perturbatively evaluate the statistics of smoothed 
field. 

• Genus curve against the scaled threshold in perturbation theory is discussed and a theory-motivated new 
asymmetry parameter in genus curve is introduced. 

This paper is structured in the following way. §2 and §3 are carrying out a mathematical exercise of expressing high- 
order terms in terms of skewness parameters; the physics of gravity only enters once one actually calculates the skewness 
parameters in §4. Thus §2 and §3 are a derivation of an extension of the Edgeworth expansion, and the basic results 
are equations (2.22) and (2.85). In §3, popular statistics of smoothed fields are examined. The second-order expressions 
in terms of the skewness parameters are given for PDF, level-crossing statistics, 2-dimensional and 3-dimensional genus 
statistics, 2-dimensional weighted cxtrema, and Minkowski functionals in §3. 1— §3.6. These quantities are re-expressed as 
functions of the volume-fraction threshold in §3.7. The skewness parameters for several cosmic fields are calculated in §4, 
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applying the second order perturbation theory. Detailed calculations of the simple hierarchical model, the 3-dimensional 
density field, the velocity field, the 2-dimensional projected density field, and the weak lensing field are given in §4.1— §4.5. 
The effect of biasing on the skewness parameters are discussed in §4.6. We discuss implications of our second-order results 
in §5. We introduce a theory-motivated asymmetry factor to characterize the shift of the genus curve in this section. 
The conclusions are given in §6. Useful Gaussian integrals including Hcrmitc polynomials are given in Appendix A. The 
bispectrum in a two-dimensional projected field is reviewed in Appendix B. The symbols in this paper are summarized 
in Appendix C. 

2. PERTURB ATIVE EXPANSION OF STATISTICS 

2.1. Smoothed fields 

We consider a cosmic random field f(x) which represents any field constructed from observable quantities of the universe, 
such as three-dimensional density field, velocity field, or two-dimensional projected density field, shear or convergence field 
of weak lensing, temperature fluctuations of CMB, etc. The coordinates x can be either three-, two-, or one-dimensions. 

We assume the field / is already smoothed by a smoothing function Wr with smoothing length R which cuts the high 
frequency fluctuations which suffer strongly nonlinear effects: 

f(x) = J d d x'W R (\x - x'\)f iaw (x'), (2.1) 

where d = 1,2,3 is the dimension of the space x, and / raw is a raw, unsmoothed field. Two of the most popular 3D 
smoothing functions are the tophat smoothing function 

W R (x) = ^Q(R-x), (2.2) 



W*W= , n J,* n * ^( (2-3) 



and the Gaussian smoothing function 

(2^) 3 / 2 i? 3 ^V 2i?V 

We also assume that the mean value of the field / is zero, and that the existence the variance <7o 2 : 

</)=0, </ 2 >=<7 2 . (2.4) 
It is convenient to introduce a normalized field a which has a unit variance as follows: 

a=—, (a 2 ) = l. (2.5) 
co 

2.2. Expressing the Non-Gaussian Statistics by Gaussian Integration 

The statistics of a smoothed cosmic field we are interested in are the functions of the field a and its spatial derivatives 
as we will see in the following sections. We denote the series of spatial derivatives by a set of variables (A^) which is 
defined by, for example, in 3D case, 

(Aft) = (a, diet, d 2 a, d 3 a, d 2 a, d 2 a, d 3 2 a, did 2 a, <9i<9 3 a, d 2 d 3 a, . . .) , (2.6) 

where di = d/dxi. For convenience, the index is denoted as fi = 0, 1, 2, 3, (11), (22), (33), (12), (13), (23), ... in such 
cases. In this example of equation (2.6), only the value of the field and of the derivatives on a single point is considered, 
but in general, more than two points can be considered. The set A^ forms multivariate random fields, which is denoted 
as an ./V-dimensional vector A in the following. The dimension N is the total number of derivatives which appear in the 
definition of the statistics we are interested in. The statistical information is described by the multivariate PDF, P(A). 
The Fourier transform of the PDF is the partition function: 

/oo 
d N AP(A) cxp(iJ • A). (2.7) 
-oo 

At this point, the cumulant expansion theorem (e.g., Ma 1985) is very useful. This theorem states that lniT is the 
generating function of the cumulants, M^"L^ n = (An ■ ■ ■ ^4 Mll ) c : 

oo . N N 

ln ^( J ) = E S E • • • E M ^n^ ■■■J.n- (2-8) 

71— 1 /il — 1 fJ^n — l 

It follows from {/) = that (A^) = 0, and first several cumulants are given by 

M« = 0, (2.9) 

M^ 2 = (A^M , (2.10) 

M £kw = (A.i4«A.s> . (2-H) 
^ii)*2M3A»4 = A^AO ~ (AnA^) (A^Am) ~ (AuiA^) (A^A^) — (AhAO (A^A^) > (2-f2) 
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and so forth. From equations (2.8) and (2.9), the partition function is given by 



^)=cxp(-ij T Mj)exp(f^ £ MiZl^-.-jA, (2.13) 
\ ' \n=3 ' m,-n n / 

(2) 

where M is an N x N matrix whose components are given by M^J . On the other hand, the equation (2.7) is inverted as 

1 r°° 

P{A) = j2^W )_J NjZ V) e *v(- iJ - A )- ( 2 - 14 ) 

Substituting J p — > id/dA^ in the last term of equation (2.13), the distribution function of equation (2.14) can be 
transformed in a form 



,«0 _ 5 " 



^) = ex p(EHf E Z.ba.. ) Pg ^ (2 - 15) 

where 



OC / 1 

■N T „ ATA 1 7T , 



P G (A) = ^-p y ^^Jexp^J.A--J 1 Mj) (2.16) 

zexpf-^M-U), (2.17) 



(2^) w / 2 VdelM V 2 

is the multivariate Gaussian distribution function characterized by the correlation matrix M. 

Any statistical quantity of a smoothed cosmic field is expressed by an average (F) of a certain function F(A) as we 
will see in the following sections. Thus, from equation (2.15), we obtain 

poo 

(F) = / d N AP(A)F(A) (2.18) 

J OO 

= (°*»[th E 8n ., A . ) F (A)\ , (2.19) 



where 



n! ^ ""^•••dA,, 



r oo 

(■■■)g= / d^P G (A)F(A) 

«/ OO 



(2.20) 



denotes the averaging by the Gaussian distribution function of equation (2.17). 

This form, equation (2.19), is useful when the deviation from the Gaussian distribution is not large. In principle, this 
equation reduces the general statistical averaging procedure to Gaussian integrations. However, it contains the infinite 
series, thus we have to truncate this expression by some criteria. In most cases of interest, the weakly nonlinear evolution 
of cosmic fields satisfy AfW ~ 0(cr o ™- 2 ), as we will see in § 4. When this relation holds, we can expand the distribution 
function to arbitrary order in er - In the following, we assume this relation, and introduce the normalized cumulants: 

-w ^ M {n) 

M (n) _ ^Mi-Mn (221) 

which are assumed to be of order one in terms of ao- In this case, the equation (2.19) is expanded as, up to 0(<7o 2 ), 
(F) = (F) G + ^ E ^™ ( F m3 ) G ^ 



-Vm (4) (f ) + - Vm (3) m (3) (f 



(To 2 + O(a 3 ), (2.22) 



where we introduce the notation, F^ lfJi2fl;i = d 3 F/dA fil dA^ 2 dA fl3 etc. The calculation of the factors (F^...)q is performed 
for individual statistic which gives the explicit form of the function F. 

2.3. Two-point Correlations 
In the expansion of equation (2.22), we need to calculate the Gaussian average of derivatives of the function F: 

(F,w~)g = J d N AP G (A)F m ..(A) = ——±j== J d^cxp (-^M-'A^j F m ..(A). (2.23) 

In most cases, the evaluation is analytically feasible because only Gaussian integration is needed. For the evaluation of 
such integration, we need the correlation matrix M. Throughout this paper, we consider statistically homogeneous, and 
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isotropic fields, in which case, the correlation matrix M is simplified because of the symmetry. In fact, the correlation 
matrix takes the following forms: 



aa) = 1, 
aa ti ) = 0, 

a ctq 



a.tOjfe) = 0, 



<?2 



(SijSki + SikSji + SuSjk) • 



(2.24) 
(2.25) 

(2.26) 

(2.27) 
(2.28) 
(2.29) 



d(d + 2) cto 2 

where ay = da/dxidxj etc. and we define the following quantities: 

ox 2 = -(/V 2 /) = -(aV 2 a)a Q 2 . (2.30) 

(T 2 2 = (V 2 /V 2 /) = (V 2 aV 2 a)(T 2 . (2.31) 

In most cases of interest, the derivatives of order higher than two do not appear in the definition of the statistics, so wc 
do not give explicit form of those correlations. 

In the following, wc use the notation rji = a t i and = a.ij following Bardeen et al. (1986). As is often the case, when 
the second order derivatives Qj appears in F as simple polynomials, the following transform is particularly useful: 

la, 2 



1 — ySija. 
da 2 



Cij — Cij 

This transform erase the correlation between a and , and the non-zero correlations are only 

(a 2 ) = 1, 

{m 2 ) = (m 2 ) = (m 2 ) = 



(Ci 2 i> = (C 2 2 2 ) - <C 3 2 3 > 



CT2 



d(d + 2) (T 2 
(C11C22) = (C11C33) = (C22C33) = 

1 (T 2 



d + 2^ 2 



1 



(T 2 



3d 

2 



d(d + 2) cto 2 
2 



d + 2 , 
—7-7 



<Cf 2 > = (CQ = (C 2 2 3 > 



d(d + 2) ao 2 ' 



where 
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(Tl 
CTQCT2 



(2.32) 

(2.33) 
(2.34) 

(2.35) 

(2.36) 

(2.37) 

(2.38) 



The Gaussian integration is straightforward if the function F is expressed by polynomials of dj . 

If the function F is more complicated in which Qj are not simply given by polynomials, it is useful to completely 
diagonalize the correlation matrix M of equations (2.24)-(2.24). We introduce the following transform, which is quite 
similar one in Bardeen et al. (1986): 



<T2 \^ (TO 2 , 



00 Cn ~ C22, 

(T 2 2 

^0 Cll + C22 — 2^33 



(2.39) 
(2.40) 

: = . (2.41) 

(72 Z 

If d = 1, wc ignore the variables y, z, and similarly, if d = 2, we ignore the variable z in the above equations and in the 
following. The above transform is similar to Bardeen et al. (1986) but notice it is not identical. This transform completely 
diagonalize the correlation matrix: 

(a 2 ) = 1, (2.42) 



(m 2 ) = (V2 2 ) = (m 2 ) = -/A 

d (To 

<.* 2 ) = l- 7 2 , 



(2.43) 
(2.44) 
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(y 2 ) 

(z 2 ) = 



1 



d(d + 2) : 



d{d + 2)' 



(C12 2 ) = (Ci3 2 > = (C 23 2 > = 



1 



ct 2 



(2.45) 
(2.46) 
(2.47) 



d{d + 2) cr 2 ' 

and all non-diagonal correlations are zero. For later convenience, we write the inverse relation of the transform of x, y, 
and z: 

(2.48) 
(2.49) 
(2.50) 



Cn = --p-{x + 4y + 2z + la ), 

4(T 
(72 

-- — (x - 4y + 2z + 7a), 

4(70 

02 



,22 



C33 = - — {x-2z + 1 a) 1 

If £33 does not appear in the function F, the variable z should be omitted in the above equations. If both £22 and (33 do 
not appear in the function F, the variables y and z should be omitted. 

After expressing the function F in terms of the diagonalized variables a, r/i, x, y, z and Qj (i < j), the calculation of 
the Gaussian integration of equation (2.23) is performed. 

2.4. Three-point Correlations 

In this paper, only first two terms in equation (2.22) arc considered. Thus we evaluate here. When the spatial 

symmetry is taken into account, the quantity reduces to the following expressions: 

(a 2 a,i) = 0, 

(a 2 V 2 a) 



d 

(a 2 V 2 a) 



a 2 a^j) = 

aa^a.j) = — 
aa^a.jk) = 0, 

. (aV 2 aV 2 a) , r . r . . . . 3((Vci'Va)V 2 a) 
aa ija fci) = , , (Oijdki + dikdji + o a o jk ) - 



2d 5ij ' 



d(d + 2) 



d(d- l)(d + 2) 



Sij Ski - 2 (SikSji + SuSjk) 



((Va • Va)V 2 a) 

Ot,iOijkCt,lm) = 0' 



SijSki — - {SikSji + SuSjk) 



(2.51) 
(2.52) 

(2.53) 
(2.54) 

(2.55) 

(2.56) 

(2.57) 

(2.58) 

(2.59) 
(2.60) 



etc. To prove the above equations, the following identities for isotropic fields are useful: 

(a^jttjj) = -^(Va • Va)V 2 a), 
3 

-(Va • Va)V 2 a) + (a(V 2 a)(V 2 a)). 

Although the above equations are valid for d 7^ 0, 1, the case d = 1 is obtained by just ignoring the terms with (d — 
Generally, more complicated quantities can be appeared for but the above relations are sufficient for our applications 

in this paper. 

According to the spatial symmetry of the above equations, the third-order correlations M^J X are explicitly given. At 
this point, it is useful to define the following quantities: 

</ 3 > (« 3 > " 



5(2) 



3(/ 2 (V 2 /)) 3 (a 2 (V 2 a))a 
4 



4 cr 2 ov 



3d ((V/-V/)(V 2 /)} 



3d ((Va- Va)(V 2 a))a a 



2(d-l) ct! 4 2(d-l) 



^1 " 



(2.61) 
(2.62) 
(2.63) 



s (2) = (/(V 2 /)(V 2 /)) = (a(V 2 a)(V 2 a))a ^ (2 ^ ] 

We call quantities skewness parameters. The first one is usually called as skewness, and the others are its 
derivatives. They are to be calculated from usual perturbation theories in § 4. Using these quantities, the third-order 
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correlations are given by 



M (3) - ^°~> 
IVI ooo — D > 

Ivl 00i — u ' 



M, 



(3) _ 4 °1 2 5 (1) 



00(ii) 

^(3) 
00(ij) 



3dcr 2 



Mi 3 L=0, (i<j), 



m(3) _ _l£L?(i) 

M$ = 0, 



(Mi), 



m' 3 ,L = 

M n (3) . u .., 



o. 



(3) 



1 ai 4 
"d 2 ((2 + 2)a 4 
1 CTl 4 



2(d- 1)S (2) -3dsf } ) , 



o(«)Ctt) d 2 (d + 2)(To 4 
= 0, (i < fc): 
M( 3 ) = 1 ^ 



2S* (2) + dS. 



(Mi), 



s< 2 > - M 2) 



(* < i) 



M, 



(3) 

o(ij)(fe;) ~ 

M. (3 L = 0, 
M (3) ,=. 
m*- 3 ** — n 



0, (i <j,k <l,i^k), 



- — — 5 (2) 
3d 2 do 4 ' 

(i < *), 



(Mi), 



M 



(3) 



ij(kk) 



M: 



(3) 



o, (t ^ i), 
i CTl 4 



(* < i), 



M 



(3) 

ij(fe0 



M 



(3) 



0, (i < j,k <l,i^k), 
0, 



2.65) 
2.66) 

2.67) 

2.68) 

2.69) 

2.70) 
2.71) 

2.72) 

2.73) 

2.74) 

2.75) 

2.76) 
2.77) 
2.78) 

2.79) 

2.80) 
2.81) 

2.82) 

2.83) 
2.84) 



i(jk)(lm) 

and so forth, where repeated indices are not summed over in the above equations, and is symmetric under permutation 
of its indices. Thus, denoting F^ v \ = (F m „a)g f° r simplicity, 

E ^S^A = M$ Q F 000 + 3 E M$ {U) F 00{U) + 3 ^ M$F 0ii + 3 £ F 0(ii)(ii) 



+ 6 E M Sii)te) J? o(«)( 3 ;) +3E M «ife)fe)^fe)fe) + 3 E M «W) F «W) +6E M «W)^W) + 



2S( 2 ) 



d cr 2 



^2 E ^bo(ii) - E ^0"^ 



01 



+ 



d 2 (d + 2)a 4 

3S< 2) ^i 4 
d(d + 2) cr 4 



3d- 



3 ( d _ : ) E F o(«)(«) ~ 6 E F o(u)(jj) + y E F °(v)(*i) 

i i<j i<j 



ijtj i<j 



3 E F 0{H){ii) + 2 E ^ ( 



0(»»)0j) +E^°( i ^ 'J' 1 



(2.85) 



This equation gives the second-order correction term of the statistical quantity (F) through equation (2.22). Once the 
field / is specified, the skewness parameters are calculated by dynamical perturbation theory of the field /. The 
remaining factors in the above equation are the Gaussian integrations of the derivatives of the function F, i.e., (F^ 1/ \)q. 
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P(f) 



5(0) 

l + ^o— H 3 (u)+O(a 2 ) 



(3.2) 



These factors can be calculated once the function F is given. In the next section, we calculate the latter factors for 
individual statistics. 

3. STATISTICS OF SMOOTHED COSMIC FIELD 

In this section we calculate the factor (F MiA , 2 ...)g for each statistic. Some of the results in this section have been 
previously presented. The Edgeworth expansion of the PDF in §3.1 is a familiar result. The result of 3D genus statistic 
in §3.4 was already given by Matsubara (1994). Wc include these old results for completeness. Other subsections present 
new results. 

3.1. Probability Distribution Function 

Perhaps the simplest yet non-trivial statistic is the PDF, P(f). The perturbative expansion of the PDF is known as 
the Edgeworth expansion (Schcrrer & Bertschinger 1991; Juszkiewicz et al. 1995; Bernardeau & Kofman 1995). As the 
simplest example, we re-derive the known Edgeworth expansion from the point of view of our general formalism above 
(see also Matsubara 1995a). 

Since the PDF is simply given by P(f) = (S(f — /))/', where 5 is the Dirac's delta function, the function F in the 
previous section for PDF P(f) is given by 

F=—S{a-f/a ) (3.1) 

Since this form of F does not depend on derivatives of a, only F oo survives in the equation (2.85). From equation (A5) 
with n = and k = 3, F 000 = {2ir) 1 / 2 e- v2 / 2 H 3 {i/), where v = f/a . Thus, the PDF is derived from the equation (2.22): 

-\/27rero 

which reproduces the well-known result. 

There is less advantage of applying our formalism to this simple statistic which can be treated by standard methods. 
Our formalism has advantages when more non-trivial statistics are considered as shown below. 

3.2. Level Crossing, Length and Area Statistics 

Next three statistics we consider here are the level-crossing statistic Ni, the length statistic iV 2 and the area statistic 
N 3 . The level-crossing statistic is defined by the mean number of intersection of a straight line and threshold contours 
of the field. The length statistic is defined by the mean length of intersection of a 2D surface and the threshold contours 
of the field. The area statistic is defined by the mean area of the contour surface in a 3D space (Rydcn 1988a; Ryden et 
al. 1989; Matsubara 1996). The level-crossing statistic is defined for ID, 2D, and 3D cosmic fields, the length statistic is 
defined for 2D and 3D cosmic fields, while the area statistic is defined only for 3D cosmic fields. For statistically isotropic 
fields, those three statistics are proportional to each other. 

In general, a statistic of the smoothed field / is a function of the threshold / t , or of the normalized threshold v = ft/cro- 
The explicit expressions of statistics Ni, N 2 , and N 3 are given by (Ryden 1988a) 

JVi(i/) = <J(a- I /)|»ji|) ) (3.3) 

N 2 { V ) = (5{a- V )[{ m f + { m f] 1/2 ), (3.4) 

NzW) = (S(a - u) [(m) 2 + ( m f + ( m f\ 1/2 ) . (3.5) 

For isotropic fields, these statistics are actually equivalent (Ryden 1988a). In fact, the distribution function of 77, = ct A 
for fixed a = v is the function of only the magnitude \f]\. Thus, using spherical coordinates for rj, one can see 

N l {v)=N 2 {v) J ^\ C os<t>\=N s {v) J ^ I sin 9 cos <j>\, (3.6) 

i.e., 

= -N 2 {v) = \N 3 {y). (3.7) 

Thus we only need to consider N\ which has the simplest expression and the rests of the statistics are automatically given 
by equation (3.7). However, if the field is anisotropic, such as the density field in rcdshift space (Matsubara 1996), the 
equation (3.7) no longer holds, and equations (3.3)-(3.5) should be used for each statistic. The equation (2.85) only holds 
for statistically isotropic fields, but the equation (2.22) is applicable even for anisotropic fields. 

Now we calculate the factor (F^ ltl2 ...)o, for the particular statistic N\. The indices fi 2 , . . . only take and 1 for N± 
statistic. Let the number of be k and the number of 1 be I. Then the factor is given by 

'' i) \ 1 ( i) " 



(F m ..) G = R(k,l) = U—) ^— ) 6{a-v)\ m \) . (:-!.<Si 
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Since the variables a and r]i are uncorrelated in the Gaussian averaging, we can use equations (A5) and (A8) in appendix A. 
Thus the above Gaussian integration results in 

l-Z 



R(k,l) 



hi-2 



cti 



Vda 



- v/2 H k {v), 



(3.9) 



where hi is given by equation (A4). Now, calculation of equation (2.85) is straightforward: 

2 



E ^S^A = SWRQ, 0) + 2S«(^J fl(i, 2 ) = i-£- ( ^ 



Vd, 



5 (0) 7I 3 (i/) + 25 (1) ifi(^) . (3.10) 



On the other hand, the Gaussian contribution is simply given by 

<^>G=i?(0,0) = i^ e - 2 / 2 . 

7T Vda 

Thus, the perturbative expansion of equation (2.22) up to second order is finally given by 

1 



n Vda 



1 



—H i {v) + —H l {v) 



<T +O{<T 2 ) 



(3.11) 



(3.12) 



The second order formulas for area and length statistics are given by equation (3.7) with the above equation. To evaluate 
the above formula, the factor and should be known. Those factors are evaluated by usual perturbation theory 
in the following section. 

3.3. 2D Genus Statistic 

The next statistics we consider are the genus statistics. The genus statistics have been attractive because it has a 
geometrical meaning of the clustering as well as the cosmological significance. The 2D genus statistic G2 is defined in 
a two-dimensional plane S in a d-dimensional space, so that d > 2 is required. In this plane S, there are contours 
corresponding to each threshold v. The 2D genus statistic is defined by the number of contours surrounding regions 
higher than the threshold value minus the number of contours surrounding regions lower than the threshold value (Adler 
1981; Coles 1988; Melott et al. 1989; Gott et al. 1990). This definition is intuitive, but an alternative, equivalent definition 
is more useful: first we set an arbitrary, fixed direction on the plane S. Then there are maxima and minima of contours 
according to that direction. These points are classified into upcrossing points and downcrossing points with respect to 
that chosen direction. The number of those points are used to define the 2D genus statistic as the following way: 

G2 = — [(# of upcrossing minima) — (# of upcrossing maxima) 

— of downcrossing minima) + (# of downcrossing maxima)] , (3.13) 
per unit area of the surface. According to this definition, the explicit expression of the 2D genus statistic is given by 

G 2 (v) = -±{5(a-v)6( m )\ m \(u). (3.14) 

For this statistic, the indices ^1,^2, • ■ • only take 0, 1, 2 and (11). Let the number of be k, 1 be l\, 2 be l 2 , and (11) 
be to. Then we need to calculate the following quantity: 




5(a - v)5(»7i)MC: 



11 



(3.15) 



Since the second derivative £n appears as a polynomial, we just use the transform of equation (2.32) . Then, from equations 
(A5), (A8) and (A9), the above equation reduces to 

1-U-li-1m 



R(k; h,h; to) = 



(2tt) 3 / 2 VvWq 



01 



-w 2 /2 



[H k+1 (u)5 m0 - H k {v)8 m i] . 



(3.16) 



Thus, from equation (2.85), 

E K\ F ** = - 2Sil) 



Vda 
<Ji 



[2R(2; 0, 0; 1) - R(l; 2, 0; 0) - R{1;0, 2; 0)] - 2S (2) 



Cl 



\/da Q 



i?(0;0;2,l) 



-u 2 /2 



(2^)3/2 VVdo-o 

On the other hand, the Gaussian contribution is simply given by 

(F} G =i?(0;0,0;0) 



S {Q) Hi(y) + AS {1) H 2 (v) + 25 (2) 



O"! 



(2tt)3/2 Wrfao, 

Thus, the perturbative expansion of equation (2.22) up to second order is finally given by 

—H i (v) + —H 2 (v) + — 



a Q +O(a 2 ) 



(3.17) 



(3.18) 



(3.19) 
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3.4. 3D Genus Statistic 



The second order formula for 3D genus statistic was already derived by Matsubara (1994), but the detailed derivation 
was omitted. For completeness, we revisit the same quantity from our general point of view here. While the 2D genus 
statistic is defined by the number of contour lines in 2D surface, the 3D genus statistic (Gott, Melott & Dickinson 1986) 
is defined by the number of contour surfaces and the number of handles in 3D space. Thus the 3D genus is defined only 
for cosmic fields of d — 3. The 3D genus statistic G3 is defined by 

G3 — [(# of handles of contours) — (# of isolated contours)] , (3.20) 

per unit volume of the 3D space. This quantity is mathematically equivalent to —1/2 times Euler characteristic of the 
contour surfaces, and thus is proportional to the total surface integral of local curvature of contours from the Gauss-Bonnet 
theorem. Although those definition is intuitive, an alternative, equivalent definition is more useful as in the 2D genus 
case. We set an arbitrary direction in the 3D space. Then there are maxima, minima and saddle points according to that 
direction. From the number of these points, the 3D genus is defined by 



G3 = — - [(# of maxima) + (# of minima) — (# of saddle points)] 



(3.21) 



per unit volume. According to this definition, the explicit expression of the 3D genus statistic is given by (Doroshkevich 
1970; Adler 1981; Bardeen et al. 1986) 



Gz{v) = -- (6(a - v)S( m )5( m )\ m \ (C11C22 - C12 2 )) • 



We need to calculate the following quantity: 

( F ,mn2---) G = R ( k : h,h, k; "in, TO22, mi 2 ) 

h 2 



d 



k 3 



d 



, da 

n £ n Fit,, 



i<j 



5(a-^)5(»?i)5(j?2)|j?3| (C11C: 



22 




(3.22) 



(3.23) 



After applying the transform of equation (2.32), the Gaussian integration of Qj fixing a is given by 

(CllC22 — Cl2' 




cri 



dotf 



H 2 (a)4 2 \{ mij }) - Fi(a)Jr({m y }) + J^({m«}) 



(2), 



r(2). 



(3.24) 



where Jm^ is defined in Table 1. Then, from equations (A5), (A8) and (A9), the equation (3.23) reduces to 



R(k; h,l 2 , h; mu,m 2 2, m 12 ) 



hi 1 hi 2 hi 3 - 2 e 



-u 2 /2 



(2^)2 \Vda a 

4 2) ({m l3 })H k+2 (v) - j[ 2 \{m ij })H k+1 {v) + J< 2) ({ mij })H k (u) 



(3.25) 



Since the 3D genus is only defined for d > 3 and our universe has the spatial dimension 3, only d = 3 is meaningful for 
actual cosmic fields. Nevertheless, we preserve the general dimension d for some flavor of generality. Thus, from equation 
(2.85), 

3 



1 



0-1 



(2tt)2 VVAto 



e - v /2 5(0)^(^ + 65(1)^3^ + 65(2)^^ 



(3.26) 



Table 1 

Definition of . For other sets of {rmj} not listed in this table, = 0. 



mn 


TO22 


TO12 


7(2) 
J 


p) 

J l 


7(2) 











1 








1 











1 








1 








1 





1 


1 











1 








2 








-2 
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The Gaussian contribution is given by (Doroshkevich 1970; Hamilton, Gott & Weinberg 1986) 

(F) G = R(0- 0,0,0; 0,0,0) = -— ^ ( e~ v " l2 U 2 {v). 
Thus, the perturbative expansion of equation (2.22) up to second order is finally given by 



G 3 {v) 



1 / <7l 



(2tt) 2 VVrf^o 



5 (o) 



—H 5 (v) + S (1) H 3 (v) + S {2) Hi(y) 



(To + O(a 2 ) 
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(3.27) 



(3.28) 



The above equation with d = 3 agrees with Matsubara (1994) 1 



3.5. 2D Weighted Extrema Density 

Next, we consider the weighted extrema density above a threshold v. The field extrema is defined to be points where 
all the first-order spatial derivatives of the field / vanish: df/dxi = 0. The weight ±1 is associated to each extrema 
according to the number of negative eigenvalues of spatial second-order derivatives of the field. When the threshold is 
high enough, the weighted extrema are approximately identified with the field peaks. 

Mathematically, this statistic is equivalent to the 2D genus statistics, according to the Morse's theorem (Morse, & Cairns 
1969; Adler 1981). Therefore, we do not need to calculate separately to obtain the result on this statistic. However, we 
present the derivation of the extrema density as an alternative calculation to the 2D genus. 

The weighted extrema in 2D field is given by 

Pc2 { v ) = (9{a - u)8( m )8(r] 2 ) (C11C22 - C12 2 )) • (3.29) 
Following the similar calculation of previous examples, and using equations in Appendix A, we obtain 

( f ,ihh2-)g = R {k; h,h; mu,m 2 2, m 12 ) 

k 2 / a \ h 2 



1) n(£) n(^) ^--)%i)^ 2 )(cnC22-c 12 2 ) 



da 



( a i 



d 



drji 



i<j 



d 



9(k 



(2tt) 3 / 2 VVdao 



h h h h e- v2 ' 2 



J i2 \{m l0 })H k+1 { V ) - J (2 \{m l0 })H k { V ) + jf^m^})^.^) 



and the perturbative expansion of equation (2.22) up to second order is finally given by 



0-1 



e- v '' 2 \H 1 (v) + 



5<°) TT , . 25^) 

HAv) H 

6 1 ' 3 



(2tt) 3 / 2 VVdao 

This result is in fact equivalent to the 2D genus statistics given by equation (3.19). 



5(2) 

H2M + -3- 



(Jo + O(a 2 ) 



(3.30) 
(3.31) 



3.6. Minkowski Functionals 

Most of the Minkowski functionals are closely related to statistics discussed above. In this subsection, we comprehen- 
sively describe the exact relation between the Minkowski functionals of a smoothed field and the statistical quantities 
considered above. 

The Minkowski functional of k = is simply the volume fraction of the excursion set K which is defined by high-density 
regions above a given threshold v. 



v^)= 1 -j K dv. 



(3.32) 



The other functionals with k — 1,2, . ..,d are defined by the integral of the curvatures on isodensity surfaces of the 
threshold v (Schmalzing & Buchert 1997; Schmalzing & Gorski 1998). In 3-dimensions, d — 3, they are evaluated by a 
surface integration averaged over whole system of volume V (Schmalzing & Buchert 1997), i.e., 



V k i3 \v) = U d 2 A{x)vf\v,x), 

V JdK 



where the local Minkowski functionals, 



(3)/ \ L 



(3)/ \ 

v 2 >{v,x) = 
v 3 '{v,x) 



6tt \R! 
1 1 



1 1 



3 '"'^ ±-kR x R 2 

1 The notations S, T, and U in Matsubara (1994) arc related to the notations here by S = S^°\ T = 25^) /3, U = S < - 2 '> /3. 



(3.33) 

(3.34) 
(3.35) 
(3.36) 
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arc defined by the principal curvatures 1/ R\ and I/R2 of the surface oriented toward lower density values. In 2-dimensions, 
d = 2, the Minkowski functionals of k = 1,2 are evaluated by a line integration averaged over whole system of 2D volume 
(surface) V (Schmalzing & Gorski 1998), i.e., 



V, 



(2) 
k 



where the local Minkowski functionals, 



H = I/ dL{x)v[ 2 \v,x), (3.37) 

V JdK 

vf\v,x) = \, (3.38) 



v¥)(v, x ) = ±-±-, (3.39) 

are defined by the principal curvature l/R\ of the line oriented toward lower density values. 

All the Minkowski functionals for a Gaussian random field are analytically derived by Tomita (1986): 

n<» - j^rwr^ (^) V'^-.M, ("o, 

where the factor ojk — 7r fe / 2 /r(fc/2 + 1) is the volume of the unit ball in k dimensions, so that ujo — 1, wi = 2, ^2 = n, 
lo 3 = 47r/3 (Schmalzing & Buchert 1997). 

It turns out that the Minkowski functionals in 2- and 3-dimensions are identical to the statistics Ni , G2 and G3 in each 
dimensions, except for normalization factors. In fact, from the Crofton's formula (Crofton 1868), the fc-th Minkowski 
functional is given by 

(d) = I dME)x (k ){KnE y (341) 

In this formula for body K in d dimensions, we consider an arbitrary fc-dimensional hypersurface E and calculate the 
Euler characteristic x^ of the intersection K n E in fc dimensions. This quantity is integrated over the space £^ of all 
conceivable hypersurfaces. The integration measure dfi k {E) is normalized to give f £ ( d ) dfik(E) = 1. From this formula, we 

can see the statistics G3, G2, N± are proportional to the Minkowski functionals of v£ d \ V^, V^ d \ In fact, x^ 3 -* is given 
by —1 times the 3D genus (or, equivalently, 1/2 times the Euler number of boundaries, x(d(K n E))), x^ is identical to 
the 2D genus, and x^ is j us t 1/2 times the number of level-crossing points (Adler 1981). Thus, Minkowski functionals 
are given by 

V W( v ) = --J*L-G 3 {v), (3.42) 
vf )(*,) = -J**—G 2 <y), (3.43) 

V± d \v) = - m{u), (3.44) 

2 Ud-iuJi 

where the boundary of the body K is identified with the isodensity contours of threshold v and the statistics on right 
hand sides are defined in <i-dimcnsions. Therefore, we have already obtained the weakly non-Gaussian expressions for 
Minkowski functionals, i.e., the Minkowski functionals of k — 1, 2, 3 arc given by the above equations and equations (3.12), 
(3.19), and (3.28). The Gaussian parts of the above equations exactly reproduce the Tomita's formula (3.40). 
The remaining Minkowski functional is the volume functional 

V^{v) = {9{y-a)). (3.45) 

In this case, from equation (A5) in Appendix A, 

R ^((^)" 9 ^) = e -^r H ^ v) - (3 - 46) 



so that equation (2.85) reduces to 



1 / v \ e-" 2 /2 5(0) 



V^ d \ V ) = -crfc \JL) + ^-H 2 (u)a + O(a 2 ). (3.47) 

The equivalent form can be obtained by integrating the Edgeworth expansion of PDF, equation (3.2). 
All the formulas of Minkowski functionals derived above for < k < 3 fit into a single expression: 

K (C V) = - — ( -^-Ve"" 2 / 2 

y k I ) (27r )(fc+i)/2 Ud _ kUJk 



l -S^H k+2 {u) + k -S^H k {u) + ^-^l S ^H k _ 2 (u) 



a Q + O(a 2 )\. (3.48) 
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3.7. Reseating the Threshold Density by Volume Fractions 

The density threshold v so far is simply defined so that isodensity surface is identified by / = vuq. However, the 
horizontal shift of the nonlinear genus curve etc. is considerably attributed to the nonlinear shift of probability distribution 
of the density field (e.g., Gott, Weinberg & Melott 1987; Matsubara & Yokoyama 1996). In order to cancel the latter 
shift, the threshold v is defined so that the volume fraction /y on the high-density side of the isodensity surface equals 
to 



fv 



i r 

v/2tt h 



dte~ t2 1 2 . 



(3.49) 



In fact, most of the work on genus analysis uses the genus curve plotted against the volume- fraction threshold v. Re- 
cently, Seto (2000) re-expressed the weakly non-Gaussian formula of genus curve (Matsubara 1994) in terms of V, using 
perturbative expansion of the probability distribution function of the density field (Juszkiewicz et al. 1995). We follow 
this method to re-express the weakly non-Gaussian formulas of various statistical quantities derived above in terms of the 
volume-fraction threshold v. 

The relation of v and v of weakly non-Gaussian field is simply given by equating the two equations (3.47) and (3.49). 
Up to first order in ctq, the relation reduces to 



5(0) 

v+—H 2 {V)<j + O{a 2 ). 
o 



(3.50) 



It is straightforward to rewrite the various analytical formulas we derived above. The results for level-crossing statistic, 
2D and 3D genus are respectively given by 

! . , - i ■ - 

3 

2 



G 2 {V) 



7T Vda 



cto + O{a 2 ) 



}■ 



(3.51) 



1 



(. 



(2tt)3/ 2 \Vdo- 
x f#i(£) + 



72 



<Jo+0{a 2 



G z iy) = - 



(2^) 2 \Vdvo 



x [H 2 {v) + - S<°>) H 3 (v) + (iS™ - 5W) Hi(V)] a a + 0(cr a 2 )} 



(3.52) 



(3.53) 



The results for Minkowski functional of k = 1,2,3 are again given by these equations and equations (3.42)-(3.44). All 
the Minkowski functional for < k < 3 fit into a single expression: 

k 



Vt d \v) = 



1 



( CT i 



(2 7 r)( fc + 1 )/ 2 w d - k u k \Vda 
k 



-572 



x \ H k _i{v) + 



(sW-SW)H k (v) 



+ 



k(k-l) 



(SW - 5 (0)) Hk _ 2[V) 



an + O(a 2 ) \ . (3.54) 



Remarkably, the highest-order Hcrmitc polynomial in the non-Gaussian correction terms vanishes in each case. In addition, 
the skewness parameters only appear as combinations of the form, — S^°\ which makes the result simpler compared 
with the original form with the direct threshold v. As we will see in the following sections, the numerical values of 
S( a \ (a — 0, 1, 2) are quite close, or even identical in some special models. This means that the non-Gaussian corrections 
of the above statistics are smaller with the rescaled threshold v than with the original threshold v. This is one of the 
central results in this paper. This tendency is in agreement with the analyses of numerical simulations. 

4. SKEWNESS PARAMETERS FOR SMOOTHED FIELDS 

We need to know the skewness parameters, for the evaluation of the second-order perturbative terms of equation 
(2.85). These quantities can be calculated by usual perturbation theory, once we specify the cosmic field, /. In each 

example of the previous section, the quantity S 2 does not appear so that we evaluate S(°\ S (1 \ and 5^ in this section. 

The other kinds of skewness parameters like can be similarly evaluated without difficulty. 

4.1. Hierarchical Model 

Before we explore actual cosmic fields, we consider a simple, phenomenological statistical model, i.e., the hierarchical 
model of higher order correlation functions (e.g., Peebles 1980). In this model, Appoint correlation function is a sum of 
N — 1 products of the two-point correlation function. Specifically, the three-point correlation function is given by 

(/(xi)/(a:2)/(aJ3)> = Q Kf(xi)f(x 2 ))(f(x 2 )f(x 3 )) + (/(x2)/(x 3 ))</(x 3 )/(xi)) 

+ </(* 3 )/(*i)></(sb 1 )/(* 2 )>]. 
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We assume the field / in the above equation is already smoothed. In this case, skewness parameters, equations (2.61)- 
(2.64) are given by straightforward calculations (Matsubara 1994): 

S<°> = = = 3Q. (4.2) 
These values depend on a hierarchical amplitude, Q, which is a free parameter of this model. The relative amplitudes 
among are not freely adjusted in the above equation. In the case of volume-fraction threshold, the first non-Gaussian 
correction of the various statistics considered in the previous section is absent, since they depend only on - . 

4.2. 3D Density Field 

Next, we consider the three-dimensional density field. The skewness parameters of this field in perturbation theory 
are already calculated by Matsubara (1994) and Matsubara & Suto (1996), using Fourier transforms of the field. We 
comprehensively review this calculation here for completeness. There is an alternative way to calculate the skewness not 
depending on Fourier transforms (Buchalter & Kamionkowski 1999). 

The cosmic field / is identified with the 3D density contrast, p/p — 1, where p is the density field. The dimension this 
field is defined in is three, d = 3. The Fourier transform of the field is useful in the following: 

f(k) = J d 3 xe- ik x f(x). (4.3) 

In this notation, two- and three-point correlations in Fourier space have the forms, 

[J(ki)f{k 2 )) = (2n) 2 S 3 (k 1 + k 2 )P(h), (4.4) 

/(fci)/(fe 2 )/(fc3)) = (2TT) 3 S 3 (k 1 +k 2 + k 3 )B(k 1 ,k 2 ,k 3 ). (4.5) 

The above forms are the consequence of the statistical homogeneity of the space, where the functions P and B are the 
power spectrum and the bispectrum, respectively. Thus, the variance and its variants [Eqs. (2.4), (2.30), and (2.31)], in 
their Fourier representation, are given by 

and the skewness parameters of equations (2.61)-(2.64) are given by 

^v/pp** 1 ^" 11 (4 ' 7) 

d 3 h d 3 k 2 
4a 2 (7i 2 J (2tt) 3 (2tt) 

* (2) - ^4 / ^^( fc i ' + k 2 fB( kl , k 2 , \ kl + k 2 \). (4.9) 

For initial random Gaussian density field, the second order perturbation theory predicts the power spectrum and the 
bispectrum as follows (e.g., Peebles 1980; Fry 1984; Bouchet et al. 1992; Bernardeau 1994a): 

P(k) = P LW (k)W 2 (kR) + OK 4 ), (4.10) 



= -4— / | fcl + k 2 \ 2 B( kl , k 2 , \ kl + fc 2 |), (4.8) 



B(ki,k 2 ,k 3 ) = 



1 + £ +( ^ + M^ + (1 _£/ fc i- fc ^ 



- • - - PLmik^P^mik^Wik^Wik^WihR) 
k\ k 2 ) k\k 2 \ k\k 2 ' 

+ cyc.(l,2,3) + O(a 6 ), (4.11) 

where Phm(k) is the linear power spectrum, and E is a weak function of cosmology (Bernardeau 1994a; Bernardeau 
et al. 1995). The field smoothing corresponds to the multiplication of the window function W(kR) which is the three- 
dimensional Fourier transform of the smoothing function Wr. It is a good approximation to use the value of E for 
an Einstein-de Sitter universe, E = 3/7, in most cases. The explicit form of the function E in terms of cosmological 
parameters Qq and Ao is given by Matsubara (1995b), which is accurately fitted by 

£Rj 3 Qo _ 1/ 3o_Ao^_3 AologioQo \ (412) 

The perturbation theory is considered to be an expansion by a parameter ao- In this respect, the power spectrum P and 
the bispectrum B is of order ao 2 and ao 4 , respectively. 

Substituting equation (4.10) into equation (4.6), we obtain 

/k 2 dk 
-^k 2 ip L1N (k)W 2 (kR), (4.13) 

up to the lowest order. Similarly, substituting equation (4.11) into equations (4.7)-(4.9), and introducing new integration 
variables, l\ = \k\\R, l 2 = l^l-R, and p = k\ ■ k 2 /(kik 2 ), we obtain 

sW(fi > " £ (5)7 (I) P - (I) <W(MS w <<..M. 



STATISTICS OF SMOOTHED COSMIC FIELDS 



15 



where a = 0, 1, 2 and 




dfi 



1 + E + 



H + (1 - E)f 



W 



h z + h z + 



VF(/!)VF(/ 2 ) 



h 2 + h 2 +hh^ a=l, 
3h 2 l 2 2 (1- ,o = 2, 



(4.15) 



So far the smoothing function is arbitrary. For a general smoothing function, the above equations can be numerically 
integrated to obtain the skewness parameters for each model of the power spectrum. For some smoothing functions, further 
analytical reductions of the above equations are possible. As a popular example, we consider the Gaussian smoothing, 
W{1) = exp(— Z 2 /2), which is frequently adopted for practical purposes. We follow the similar technique of Lokas et al. 
(1995), in which they derived the skewness of the density field with Gaussian smoothing. 
For the Gaussian smoothing, 



W (Vh 2 + h 2 + 2hl 2 n) 



-Z1Z2M 



W{h)W(h) 

In this case, the following formula of the modified Bcssel function I v (z) is useful: 



/_' 



dnPi{n)e-»* = -I l+1/2 {z) 



(4.16) 



(4.17) 



where Pi is the Z-th Legendre polynomial. From this formula, the angular integration of ^ in equation (4.15) can be 
analytically performed and the result is 



5(0) 




3 fh I 



(2 + E)I 1/2 (hl 2 ) - - -1 + -i I 3/2 (hl 2 ) + (1 - E)I 5/2 (hl 2 ) 



2E 



2 \h I 

2 



3(9 + E) l 2 W 



10 



h/ 2 {hh) 



+ 



E 



!■> 2C.1V2) - ' - in --f7/2Gl^2)| 



= V2^(hl 2 ) 3 / 2 



3(3 + 2^) 9 
z h/ 2 {hh) - -7; 



3(3 + AE) 
7 



h/ 2 (hl 



10 \h 
9 

To 



h h 
'2 



10 

^3/2(^1^2) 



h/ 2 {hl 2 ) - 



18(1 — E) 
35 



-^9/2(^2) 



(4.18) 



(4.19) 



(4.20) 



At this point, it is useful to define the following quantity: 



V2~7T ( <Jq 



a+/3-2 



h 2 dk l 2 2 dl 2 

27T 2 i? 3 27T 2 i? 3 



h 



Ptjn t Pi 



R 



LIN 



|) e-^-^l^l^I^ihh). (4.21) 



In the above equation, variance parameters ao, ci are given by equation (4.13). The nonlinear correction for Oj is not 
needed because our estimate of is only lowest order in a . When the higher order corrections, e.g., third-order 
perturbation corrections, are estimated, one should be sure that all the necessary nonlinear corrections are properly taken 
into account. With this quantity, the skewness parameters are given by 



S (0) (R) = (2 + E)S^ - 35? 2 + (1 - E)S\\ 



S {1) (R) = 
S<V{R) = 



5 + 2E 



3 + 2E 
15 



o33 
D 



9 + E 



- q24 



sf - + 2(2 - E) s? 



E 



,22 



3 + 4B, 
21 1 



<33 , ± q24 
•2 + g^3 



2(1 ~ E) 
35 



c33 



(4.22) 
(4.23) 

(4.24) 



Thus the lowest order estimates of skewness parameters are given by the above equations (4.21)~(4.24). For each given 
power spectrum, the integration of equation (4.21) is straightforward. 

The resulting skewness parameters are independent on the amplitude of the power spectrum. When the power spectrum 
is given by a CDM-like model, P\,m{k) oc kTcnM 2 (k/T), where 



Tcdm(p) = In(1 2 + 3 4 P 34P) [1 + 3.89p+ (16.1 P ) 2 + (5.46p) 3 + (ejlp) 4 ]" 1 



/4 



(4.25) 



is the CDM-like transfer function fitted by Bardeen et al. (1986), and T is the shape parameter of this model, then the 
skewness parameters are functions of TR. In Table 2, we give the values of skewness parameters for CDM-like models 
for several values of TR. In this table, the value of E is approximated by 3/7. Since the skewness parameters are weak 
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functions of TR as seen from the table, one can interpolate the values in this table to obtain the values of arbitrary scales 
for practical purposes. 

When the power spectrum is given by a power-law form, 

P Lm (k)=Ak n °, (4.26) 

the integration of equation (4.21) can be analytically performed. First, the simple Gaussian integration gives 

3 ~ —0, D „-4-2,-4-3 r 1 O ) ' i 4 ' 27 ) 



Second, we expand the modified Bessel function as 



r=0 



Then the equation (4.21) reduces to 

QaP _ 2 fn s + 3\ 1 - {a+fi)/2 fn s + 3\ (n s + 3 



(2m+ 1)!! V 2 



(a+m-l)/2 \ * / (/3+m-l)/2 



/n s + a + m + 2 n s + (3 + m + 2 3 1\ 

x F ^ , , m + -; - j , (4.29) 

where (a)„ = T(a + n)/Y(a) = a(a + 1) • • • (a + n — 1), and F is the Gauss hypergeometric function: 

n«,A7;^) = E^f k S- (4 - 30) 

The equations (4.22)-(4.24) and equation (4.29) give the skewness parameters S^ a \ The following recursion relations for 
hypergeometric function 

aF(a + 1, p, 7 + 1; z) = jF(a, (i, 7; z) + (a - 7 )F(a, /?, 7 + 1; 0), (4.31) 
a/3F(a + 1, + 1, 7 + 1; z) = 2llzi2 [p(a, p, 7 - 1; z) - F{a, (3, 7; z)\ , (4.32) 

(7 - a)(7 - /3)*F(a, /3, 7 + 1; z) + 7 [(a + /3 - 2 7 + l)z + 7 - l]F(a, (3, 7; «) 

+7(7 - 1)(* - /3, 7 - 1; *) = 0. (4.33) 

simplify the result: 

c(0)r \ c(Dr \ OI ? f n s + 3 n s + 3 3 l\ f n s + 3 n s + 3 5 l\ 

= >(n s ) = 3F I —j—, — J— . 2' 4 ) ~ ( ™ s + ^ ' ~~ 2~ ' 2' 4 J ' ( ^ 

0(27 \ oi7^ n s + 5 n s + 5 5 1\ 3 + 5 n s + 5 7 1\ 

5 } K) = 3F (-J", -a". 2' ij - 5 K + 4 - 4S)F (-2-' -2~' 2 5 4 J • (4 ' 35) 

In this power-law case, skewness parameters do not depend on scales R but only on power-law index, n s . Incidentally, 
S^(n s ) and S^(n s ) are identical. This is just the coincidence and is not generally the case when the power spectrum is 
not given by power-law. Several numerical values are shown in Table 3, where E = 3/7 is assumed. 

4.3. 3D Velocity Field 

Next, we consider the 3D velocity field as the cosmic field /. Since the rotational components of the velocity field 
are decaying modes of gravitational evolution in perturbation theory, we only consider the rotation- free component. The 
cosmic field / is identified with the dimensionless scalar field, 

f(x) = lv • v(x), (4.36) 



Table 2 

3d skewness parameters of gaussian-smoothed density field for cdm-like models, which are the functions of 
the product of shape parameter v and smoothing length r. the parameter e is set as e = 3/7. 



TR 


1.0 


2.0 


4.0 


8.0 


16.0 


32.0 




3.678 


3.500 


3.332 


3.201 


3.115 


3.065 




3.757 


3.566 


3.377 


3.228 


3.129 


3.072 




3.657 


3.662 


3.701 


3.783 


3.903 


4.046 
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where H is the Hubble parameter. One can also consider other quantities like radial component of the velocity field, 
V = n ■ v, where n is the line-of-sight normal vector. Those quantities are more complicated than a simple divergence. 
We illustrate only the simplest case in this paper. The second order perturbation theory predicts power spectrum and 
bispectrum of the velocity field as follows (e.g., Bernardeau 1994a): 

P(k) = g f 2 P L1N (k)W 2 (kR) + O(a 4 ), (4.37) 

2" 



B(k u k 2 ,k 3 ) = -g f 3 



1+E„ 



fci 



ki\ fci • fc 2 



k 2 J kik 2 



+ (1 - E v ) 



ki ■ fci 

k\k 2 



+ C yc.(l,2,3) + O( ( 7o 6 ). 
In the above equation, the factor gf is a logarithmic derivative of the growth factor, given by 

5/("o, A ) 



J PLIN(fcl) J PLIN(fc 2 )IF(fcii?)M/(fc 2 i?)W(A:3i?) 

(4.38) 



din a 



70 V 2 



(4.39) 



(Lightman & Schcchtcr 1990; Lahav et al. 1991) and D is the linear growth factor and a is the expansion factor. The 
logarithmic derivative is evaluated at present. The factor E v is a weak function of cosmology. It is a good approximation 
to use the value for an Einstein-de Sitter universe, E v = —1/7. The explicit form of the function E v in terms of fi and A 
is given by Matsubara (1995b), which is accurately fitted by 



El » + 1 ^ 3 -11/200 A o (, 7 \ 
—2~ -_^1--A log 10 f7 j. 



(4.40) 



The similarity of the equations (4.37) and (4.38) for velocity field to the equations (4.10), (4.11) for density field is 
obvious. We can easily see the skewness parameters for the velocity field are obtained by similar equations as (4.22)- 
(4.24): 

S<®(R) = -— [(2 + E^S' 1 - 3S , 2 + (1 - E V )S?] , 



5 (1) (i?) 



9f 
3 



b + 2E v 



cl3 



9 + E v 



S f - + 



2(2 — Ey) 13 1 - E v 



Q 22 

°3 



S<»{R) = - — 
9f 



j~ 2Ey „ 33 1 „ 2 4 3 + AE V 33 1 24 2(1 E v ) 3 3 

15 5 bl ~^l~ b2 + 5 ' 3 35 bi 



(4.41) 
(4.42) 
(4.43) 



where S^f is given by equation (4.21) without any modification. In Table 4, we show the values of skewness parameters 
of the velocity field for the CDM-like models for several values of TR. In this table, the value of E v is set as —1/7. 
For the power spectrum of the power-law form, 



S^(n s )=S^(n s 

S( 2 Hn s ) = -- I 
9f 



1 

9f 



3F 

5 n s 



n s + 3 n s + 3 3 1 
2~ ' 2~ ' 2' 4 
5 5 1 
~ 2 ! 4 



-{n s + A-AE v )F 



- (n s + 2 - 2E V )F 
; + 5 n 



n s + 3 n s + 3 5 1 
2~ ' 2~ ' 2' 4 
+ 5 7 1 
2~ ' 2' 4 



(4.44) 
(4.45) 



Several numerical values are shown in Table 5, where E v = — 1/7 is assumed. 

4.4. 2D Projected Density Field 

The projection of the density field p p defines the 2D cosmic fields on the sky (d = 2). Here, we derive the skewness 
parameters for this field. The 2D projected density field with top-hat kernel is investigated by Bernardeau (1995). Here, 
we are interested in Gaussian kernel for our purpose. 

In a Friedman-Lemaitre universe, the comoving angular diameter distance at a comoving distance \ is given by 

' sinh (xV^K) 



Sk(x) 



X, 
sin 



(xVk) 



-,K<0, 
K = 0, 

K>0, 



(4.46) 



Table 3 

3d skewness parameters of gaussian-smoothed density field for power-law models, which are the functions of 









THE 


SPECTRAL 


INDEX n. 








n 


-3.0 


-2.5 


-2.0 


-1.5 


-1.0 -0.5 


0.0 


0.5 


1.0 


5(2) 


4.857 
3.815 


4.400 
3.720 


4.022 
3.665 


3.714 
3.652 


3.468 3.280 
3.680 3.754 


3.144 
3.877 


3.061 
4.054 


3.029 
4.294 
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depending on the sign of the spatial curvature K = i? 2 (£!o + Ao — 1). Thus, in spherical coordinates, projected density 
field in 2D is given by 

Pp (9, <j>)= J d X S K 2 (x)n( X )p( X , 0, <p; r - X ), (4-47) 

where n(x) is the selection function without volume factor, normalized as / dxS^(x) n ix) = 1) an d p(x,S,4>] T ) is the 3D 
comoving density field 2 . The present value of the conformal time r — j dt/a is To- 

The projected density contrast is defined by p p /p p — 1, where one can see p p = p. We identify the projected density 
contrast with the 2D (d = 2) field /. Since the smoothing angle Of is much smaller than tt in most of the interested cases, 
we consider the small patch of the sky of the vicinity of the polar axis, 9 <C 1. With this approximation, we introduce the 
variables 9i = 9cos<j), and 9 2 = #sin</>, which are considered as 2D Euclidean coordinates, 6. Therefore, the projection 
equation is given by 

f(6) = j d X S*(x)n( X )5Mx, Sk(x)0; t - X ), (4.48) 

where 630 (x,t) = p/p— 1 is the density contrast at comoving coordinates x and conformal time r. 

The power spectrum and the bispectrum for the above projected field are given by the Limber's equations (B2) and 
(BIO) in Appendix B, with q(x) — n (x) : 

^2dM = j d X Sk(x)n 2 (x)P3B(^^ y To - X), (4-49) 

/2 3 w l ^2 ^2 

dxS K {x)n i {x)Bm{ Sk ^ , g^y S K (x) ' T0 ~ X ^ 

where P, B are 2D projected power spectrum and bispectrum, respectively, of the field /, and P 3 d, -B3D are 3D power 
spectrum and bispectrum, respectively. The 3D power spectrum and the 3D bispectrum are evaluated by the second order 
perturbation theory. They are similar to equations (4.10) and (4.11), but we have to take into account the time-dependence 
here. They are given by 

P3D(k;T -x) = D 2 (x)PuN(k) (4.51) 

2" 



B 3D (ki, k 2 , h; tq-x)=D (x) 



i + E( x ) + 7- + tt -in— + (! - E (x)) 



PLIN(fcl)P L IN(fc 2 ) 



ki k 2 J kik 2 \ k\k 2 

+ cyc.(l,2,3), (4.52) 

where D(x) is the linear growth factor at conformal lookback time (i.e., at conformal time r = r — x), which is 
normalized as D(Q) = 1. The following fitting formula (Carroll, Press & Turner 1992) is useful: 

an ^o 4/7 -A + (l + »o/2)(l + A /70) 

~ fl ' f7 4 / 7 -A+(l + fi/2)(l + A/70) ' { ' ' 

where f2 and A are time-dependent cosmological parameters at conformal lookback time X- The variable E(x) is a weak 
function of time and cosmology, and for Einstein-de Sitter universe, E = 3/7. This quantity E is the same we used in 3D 
density field, but here we also take into account the time-dependence. It is accurately approximated by 

E » 'in- 1 / 30 - A ( i _ hio gw n) . (4.54) 



7 80 V 2 

The variance parameters of the smoohted projected field are given by 



where 



= j ^-^P^W 2 ^) = J d X S 2 (x)n 2 (x)D 2 (x)^ 2 [S K {x)0f] , (4.55) 

E/(i?) = R 2 ^ 2 J ^k 2 ^P LW (k)W 2 (kR) = J ^/ 2j >lin W 2 (l). (4.56) 



2 The comoving density field is defined by the density per unit comoving volume and thus satisfy p = constant. 



Table 4 

3d skewness parameters of gaussian-smoothed velocity field for cdm-like models. the factor —qf is 
multiplied and the values are almost independent on cosmological parameters. 



TR 


1.0 


2.0 


4.0 


8.0 


16.0 


32.0 


-9fS<»> 
-9fSV 


2.456 
2.551 
1.963 


2.232 
2.319 
1.877 


1.995 
2.065 
1.779 


1.773 
1.826 
1.683 


1.581 
1.620 
1.601 


1.420 
1.448 
1.530 
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The skewness parameters of the smoothed projected held are given by 
S (a) (6 f ) 



where 



and 



ao^f 4 \c7i9f 

cW < r »^(It)7 



J d X S K 2 ( X )n 3 ( X )D\x)Zo A - 2a [S K (x)0f] Si 2a [S K (x)0f] C (a) [S K { X )9 { 



lidli l 2 dl 2 
~~2^ 2^T 



1 + E- 



PlAN 4 PT 



R 



h h 



LIN 



l -l)w 2 (h)w 2 (i,u (/,../,). 



(4.57) 
(4.58) 



/x + (1 - E)f 



W{h)W(l 2 ) 



V 7 ! - I*' 
2, a = 0, 

x { h 2 + l 2 2 + hhn, o=l, (4.59) 
4/ 1 2 ? 2 2 (l- A1 2 ),a = 2, 

So far the smoothing function is arbitrary. For a general smoothing function, the above equations can be numerically 
integrated to obtain the skewness parameters for each model of the power spectrum. In the following, we adopt the 
Gaussian smoothing function, W{1) = cxp(— l 2 /2). For this smoothing function, the equation (4.16) holds even for this 
2D case. In this case, the following integral representation of the modified Bessel function I„(z) for v = is useful: 



h{z) = 



1 



V^ 2 



-ZfJ, 



Actually, the derivatives of the above equation, 



d/i 



(4.60) 



(4.61) 



are sufficient to perform the angular integration in equation (4.59). Moreover, one can use the property of the Bessel 
function, I' a = I\, and I' m = (I m -i + I m+ \)/2 to obtain formulas, 



1 f 1 dn 



fj,e 



1 



, - 7 ^=^ 2 e-^^ 1 -I Q (z)+ 1 -I 2 (z\ 
n J-i \fl- H 2 2 2 



i ^y.-*-|/ lW -ij iW , 



i r 1 dv 
* J-i JT^ 



From these formulas, equation (4.59) reduce to 

3 



cm 



h h 
h h 



{3 + E)I (hl 2 ) - 2 
{± + E)(l 2 + l 2 )h(hl 2 )- 



h(hh) + {l-E)h{hh 

3 



h h 



h{hh) 



+ (2 - E)(l 2 + l 2 )I 2 {hl 2 ) - ^-^hkhihk) 



= - [(5 + 3E)l 2 l 2 I Q (hl 2 ) - 2hl 2 (l? + h 2 )h(hh) - 4(1 + E)l 2 l 2 I 2 (hl 2 ) 
+ 2l x l 2 {ll + li)h{hh) - (1 - E)l 2 l 2 h{hl 2 )} . 



(4.62) 
(4.63) 
(4.64) 
(4.65) 

(4.66) 

(4.67) 
(4.68) 



Table 5 

3d skewness parameters of gaussian-smoothed velocity field for power-law models, which are the functions 

OF THE SPECTRAL INDEX fl. 



n 


-3.0 


-2.5 


-2.0 


-1.5 


-1.0 


-0.5 


0.0 


0.5 


1.0 


-9fSW 


3.714 
2.333 


3.250 
2.170 


2.848 
2.026 


2.498 
1.900 


2.191 
1.788 


1.918 
1.687 


1.670 
1.595 


1.441 
1.509 


1.222 
1.425 
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At this point, defining 

the equation (4.58) reduces to 



rhdh 

J 2^ 



_/-_/- 



h a 1 h 13 1 Im(hh), 



2(4 + E)C'q - 



4Cf 
15 



{l-E)C l 2 1 } 



E 



E 



C {2) - \ [(5 



3E)C 3 3 



2 

4Cf 



4(1 + E)Cf + 4Cf 4 - (1 - £)Cf] . 
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(4.69) 

(4.70) 
(4.71) 
(4.72) 



The two-dimensional integration of equation (4.69) is performed only once as a function of R. The result is stored as a 
table, and is used in one-dimensional numerical integration of equation (4.57) for finally obtaining the skewness parameters 
in 2D projected density fields. Functions for CDM-like models are given in Table 6 which are functions of TR. 

When the 3D power spectrum is given by a power-law form of equation (4.26), the integration by wave length l\, l 2 
can be analytically performed as in the 3D case. In fact, the parameter Ej of equation (4.56) with Gaussian smoothing 



W(l) 



-l 2 /2 



is given by 



E/(i?) 



A 



and the variance parameters are given by 



n s + 2j + 2 



47r f n s +2 j+ 2 

With similar technique used in 3D density field, 



AttR"- 
n s + 2j + 2 



)/ 



d X S K 2 - n °(x)n 2 ( X )D 2 (x)- 



(4.73) 



(4.74) 



C af3 = 



1 



2 m m! 



n + 2 



xF 



l-(a+/3)/2 



n + a - 



n + 2 



1 n 



(ct+rr 



-l)/2 

m 4 



1 



(/3+m-l)/2 

V 



2 1 2 

The equations (4.70)-(4.72) and the above equation finally give the values of 
(4.31) and (4.32) simplify the result: 



C(°'(n s )=C«(n s ) = 3F 



2 n s + 2 1 



3 



.^T )-^(n s + l-E)F 



(4-75) 

The recursion relations of equations 
n„ + 2 n s 



2 1 
■ ,2; - 

2 ' ' 4 



C( 2 >(n s ) = 3F 



n s + 4 n s +4 2 ^ 1 
2 ' 2 ' ' 4 



3(n s + 3 - ZE) 



F 



n s + 4 n s + 4 1 
2 ' 2 ' '4 



(4.76) 



(4.77) 



These functions are independent on scale R in the power-law case, but dependent on power-law index, n s . Again, C^(n s ) 
and C^(n s ) are identical only in the power-law case. Numerical values are given in Table 7, where E = 3/7 is assumed. 
For the power-law case, the skewness parameters of equation (4.57) reduces to the following simple form: 

/ d X [S K { X )f- 2n ° n 3 (x)D i ( X ) 

s (a) K) = -rz — c (a) K). 



dx[S K {x)Y~ ns n 2 ( X )D 2 (x) 



(4.78) 



It is interesting to compare the results with those of top-hat smoothing. According to Bcrnardeau (1995), the top-hat 
smoothing gives C^(n s ) = 36/7 + 3(n s + 2)/2 for power-law case. Comparing this expression with our Table 7, they 
roughly agree with each other, and have similar behavior with spectral index. In detail, the Gaussian smoothing gives 
more or less larger values. For example, top-hat smoothing gives C(°) = 5.14,3.64,2.14 for n s = -2,-1,0, respectively, 
while Gaussian smoothing gives C^ ' = 5.14,3.89,3.01, respectively. It is reasonable that the of top-hat smoothing 
differs from that of Gaussian smoothing for large spectral index, because the top-hat smoothing gathers more power on 
small scales than Gaussian smoothing. 



TR 


0.25 


0.5 


1.0 


2.0 


4.0 


8.0 


16.0 


32.0 


64.0 




4.544 
4.885 
3.723 


4.311 
4.614 
3.660 


4.020 
4.279 
3.589 


3.685 
3.890 
3.518 


3.345 
3.492 
3.467 


3.044 
3.143 
3.453 


2.797 
2.863 
3.479 


2.613 
2.653 
3.532 


2.504 
2.523 
3.584 



Table 6 

Functions C (a) (i?) for CDM-like models, which are the functions of the product of shape parameter r and 

SMOOTHING LENGTH R. THE PARAMETER E IS SET AS E = 3/7. 
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4.5. Weak Lensing Field 

The local convergence field of the weak lensing is commonly used for studying the large-scale structure of the universe 
(e.g., Kaiser 1998; Bartelmann & Schneider 1999). Assuming the situation where Limber's equation (see Appendix B) 
and also the Born approximation (Kaiser 1998) hold, the following correspondence between the convergence field k and 
the 3D density contrast S^b is useful (e.g., Mellier 1999): 

k(0) = '±H 2 n r dx'S K \ X ')n{x') [ X dx SK{ * )S * {} f~ X) 6 3B [ X ', 0S K (x'y,T X '} , (4.79) 
2 Jo Jo a{x)S K {x) 

where a( X ) is the scale factor at conformal lookback time X = r — r. The above equation is reduced to exactly the same 

form as the projected field of equation (4.48), but with the substitution, 

»<*> - -<*> - i^r r <48o) 

However, one should note that this substitution is valid only under the assumption of the Born approximation. Although 
the effect of the Born approximation on the skewness is known to be weak (Bcrnardeau et al. 1997), the validity of the 
Born approximation in general situation has not been tested in detail (Mellier 1999). There is some subtlety which could 
arise when various combinations of skewness are considered. 

4.6. The Biases 

The above expressions of the skewness parameters are for unbiased fields. The skewness parameters for biased fields are 
non-trivial. They depend on the details of the biasing scheme in the real universe which is poorly known so far. However, 
the skewness parameters are well-defined quantities, so that they are calculated from the first principle once the biasing 
scheme is given. 

Perhaps, one of the simplest, yet non-trivial case is the local, deterministic biasing. In this case, we can follow Fry & 
Gaztanaga (1993) to obtain the perturbative expansion of the biasing: 

f g = 6«+£(^ -<«»»+..., (4.81) 

where 5 g and S are the galaxy and mass density contrast, respectively. The spatial dimension is arbitrary, so that 5 g , 5 
can be functions in either ID, 2D, or 3D space. The biased variance in lowest order is given by 

<7 , g - bcr Q . (4.82) 

It is also straightforward to calculate skewness parameters. After some algebra, all skewness parameters are shown to 
transform in a same way: 

S^ = S ~^ + % (4.83) 

irrespective of spatial dimensions, d, and of kinds of skewness parameters, a = 0, 1, 2. In this framework, the parameters 
b and bi are needed. A possible way to determine these parameters is to measure the variance cr g 2 and the skewness S g 
from the observation. Theoretical models predict the variance cr m 2 and the skewness S m of the mass distribution. We 
obtain biasing parameters by b = a g /a m and 6 2 = b(bS g — 5 m )/3. The derivatives of the skewness, Sg~\ S g 2 \ are then 
obtained from equation (4.83). 

In case of level-crossing, and genus statistics as functions of scaled threshold u, the second-order corrections in equation 
(3.51)-(3.53) only depends on the difference of the skewness parameters times the variace. Quite remarkably, this type of 
second-order correction term does not depend on bias parameter at all: 

4 a) - S ( g A <7 0;g - [S (a) - S (0) l an (4.84) 

Thus, the second-order nonlinear corrections for level-crossing and genus statistics of locally biased field are exactly the 
same as that of unbiased mass density field. This can be considered as an advantage of the scaled threshold v in these 
statistics, because the biasing is one of the most embarrasing uncertainty in the analysis of galaxy distribution. However, 
one should note this result is derived under the local biasing scheme. So far the locality of the bias is not guaranteed in 
general. 





-3.0 


-2.5 


-2.0 


-1.5 


-1.0 


-0.5 


0.0 


0.5 


1.0 




6.949 
4.090 


5.965 
3.863 


5.143 
3.687 


4.457 
3.560 


3.885 
3.478 


3.409 
3.443 


3.014 
3.445 


2.688 
3.518 


2.421 
3.635 



Table 7 

Functions C (a) for power-law models, which are the functions of the spectral index n. The parameter E is set 

as £ = 3/7. 
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Stochastic argument (Dekel & Lahav 1999) of the biasing is more complicated for derivative skewness and than 
for usual skewness S(°\ because they involve the correlation between field derivatives. The phenomenological nature of the 
stochastic biasing requires many parameters which is not calculable from first principles. Therefore, the stochastic biasing 
shcme does not effectively work for our problem. More physical treatment of the biasing schemes of galaxy formation is 
needed, in which case the nonlocality of the bias could also be important (Matsubara 1999). 

5. IMPLICATIONS OF SECOND-ORDER RESULTS 

The perturbative calculations offer valuable aspects of weakly non-linear evolution of the various statistics without 
laborious parameter survey by numerical simulations. In this paper, we obtain the lowest nonlinear corrections to relatively 
popular statistics of smoothed cosmic fields. Using these results, it is interesting to see how the weakly nonlinear effect 
tends to distort the Gaussian prediction for those exemplified statistics. 

In Figure 1, various statistics for 3D density field is shown. The amplitude of each statistic is appropriately normalized 
as we are interested in the deviation from the Gaussian prediction. If we neglect the normalization, Minkowski functional 
of k = 1, 2, 3 are equivalent to the statistics Ni, G2, and G3, respectively, so that they degenerate in this figure (we should 

note the sign of V3 is inverted) . 

The rms a , which is considered as a weakly nonlinear parameter, is set ct = 0.3. A limit er — ► corresponds to the 
prediction of the linear theory, which is given by thin solid lines in the figure. This linear prediction is equivalent to the 
Gaussian fluctuations, because we assume the initial density field is random Gaussian. 

In general, the curves of statistics plotted against the direct density threshold, v (dotted lines), exhibit considerable 
deviations from Gaussian predictions. The overall tendency does not depend much on the shape of the spectrum we 
consider here, i.e., power-law spectrum with index —2, —1, 0, and CDM model with smoothing length R = 4/T, where T 
is the shape parameter of the CDM spectrum. They are consistent with the so-called meat-ball shift, which means that 
there are more isolated regions in a nonlinear field than in a Gaussian field for a fixed threshold. In fact, N±, G2, — G3 
of high value of threshold, e.g., v ~ 2, virtually correspond to the number of isolated regions, and each figure shows that 
the number is indeed increased by weakly nonlinear evolution. 

The weakly nonlinear formula for the genus curve against the density threshold, Gs(v), which was first derived by 



T 




-2 2 -2 2 -2 2 -2 2 
v, v v, v v, V v, v 



Fig. 1. — The 3D genus, G3, the 2D genus G2, the level-crossing statistic N\ and the Minkowski functionals V± of the 3D density field. All 
curves are appropriately normalized. The variance is set as ctq = 0.3. Solid lines: Gaussian predictions, dotted lines: second order predictions 
in terms of density threshold, v, dashed lines: second order predictions in terms of volume- fraction threshold, v. The initial density fluctuation 
spectrum is given by the power-law with n = —2, —1,0 and also by the CDM-like model with smoothing length R = 4/T (see text). 
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Matsubara (1994) has been compared with numerical simulations in literatures. Matsubara & Suto (1996) shows the good 
agreement of the analytic prediction with the simulations for various spectra. Colley et al. (2000) compared the prediction 
of the genus curve against direct v with the simulated SDSS data. Unfortunately, in their published paper, they have 
transcribing errors, which incorrectly made the perturbation theory considerably disagree with their data. There are also 
ambiguity in their comparison on biasing which alter the values of skewness parameters. One can guess the biasing effect 
on skewness parameters by equation (4.81). They chose the peak particles as galaxies, and the linear biasing parameter 
is inferred as b <~ 1.3, but the nonlinear parameter b 2 is not obvious in their work. Some literatures indicate the skewness 
of peaks are roughly given by 1-2 (Watanabe, Matsubara & Suto 1994; Plionis & Valdarnini 1995) but for highly 
biased peaks 6^2. If we adopt 62 = —0.5, the skewness is given by — 1.8 which is not unreasonable. If it is the 
case for their simulation, the perturbative prediction and their data completely agree with each other. The x 2_v alue per 
degrees of freedom reduces to only 1.03 (private communication with W. N. Colley & D. H. Weinberg). Obviously, we 
have to further investigate the biasing effects in numerical simulations to obtain a conclusive result. 

Most of the topological analyses of the previous work use the scaled threshold v. The dashed lines in Fig. 1 shows the 
corresponding curves. The deviations from the linear theory is dramatically reduced. This fact is empirically known by the 
analyses of numerical simulations (Gott, Melott & Dickinson 1986; Gott, Weinberg & Melott 1987). The reason for this 
reduction is mathematically due to the closeness of the values of skewness parameters S^ a \ (a = 0, 1,2), since all the terms 
of the nonlinear corrections in equations (3.51)-(3.53) depend only on — (a = 1,2). For the hierarchical model of 
equation (4.1), they are exactly zero, that means there is not any (second order) nonlinear correction for the hierarchical 
model. Since the hierarchical model is known to roughly approximate the nonlinear evolution, it is not surprising that 
more realistic fields have only small corrections of nonlinearity if they are plotted against the volume-fraction threshold, 
v. 

For power- law models, and are exactly the same. That makes the nonlinear correction for N\ or exactly 
vanishes. Thus, the nonlinear corrections for other statistics are arisen by the difference between and S^ 2 \ For the 
CDM model, there still is a difference between and but it is relatively small as seen from the Table 2. 

As for the topological statistics, G2, G3, V^ 3 ' and v£ 3 \ deviations of curves against v from the linear theory prediction 
depend on the underlying spectrum through differences of skewness parameters. The weakly nonlinear effect for redder 
spectrum of n s = —2 induces a sponge-like shift, that means the number of holes in isolated regions increases. On the 
other hand, the bluer spectrum of n s = indicates a meat-ball shift. These tendencies are qualitatively in agreement with 
the numerical results (e.g. Ryden et al. 1989; Melott et al. 1989; Park & Gott 1991). 

It is not trivial to estimate errors expected in observationally estimating statistics of smoothed cosmic fields. The 
observational errors mainly consists of the shot noise and the cosmic variace. The systematic comparison with simulations 
should be used for the error estimates. Unfortunately, previous earlier simulations does not have enough resolution to 
be quantitatively compared in weakly nonlinear regime. Canavezes et al. (1998) use much larger simulations than those 
earlier work, and gives the genus curve in weakly nonlinear regime in their Fig. 9. Relative depths of left and right 
troughs in the genus curve show slightly meat-ball shift, which is consistent with our prediction. Although their variation 
of smoothing length is rather limited, their plots ensures the slight meat-ball shift predicted by the perturbation theory 
can definitely be observed. Colley et al. (2000) plot the 2D genus curve, but they use the thin slice which suffers strongly 
nonlinear effects so that the quantitative comparison is not appropriate with weakly nonlinear results. They report the 
meat-ball shift in 2D genus, which is the same direction of the weakly nonlinear correction besides amplitude. We find 
there is stll neccesity of proper, systematic comparison between the perturbation theory and numerical simulations in 
weakly nonlinear regime. 

The decisive factor between sponge-like shift and meat-ball shift is the sign and amplitude of — S^°\ since 

5(1) _5(0) 

is exactly or approximately zero for wide range of models like power-law, or CDM-likc models. If the factor — is 
positive, the meat-ball shift takes place. If this factor is negative, the sponge-like shift occurs. For power-law case, that 
factor is positive for n s > —1.4, and is negative for n s < —1.4. The amplitude of this factor times the amplitude of the 
nonlinearity parameter, (S^ — S^)ao determines the amplitude of the meat-ball shift in the genus curve. 

Prominent shifts are seen around the two troughs of the genus curve oXv = ±v/3. At these troughs, the factor -S^ 
does not contribute to the genus curve since the accompanying factor H^{y) = v 3 — v vanishes. Therefore, one can almost 
completely characterize the shifts by a factor (S^ 2 ' — S^°')aa 7 which we call the 'genus asymmetry parameter', around the 
troughs. We propose this genus asymmetry parameter as a theory-motivated parameter of the genus asymmetry. This 
factor is observationally determined by fitting the shape of the genus curve by a form, 

G{v) oc H 2 {v) + AHz(V) + BH^v), (5.1) 

where B is the genus asymmetry parameter. In this fitting, A is much smaller if the non-Gaussian features are from purely 
gravitational evolution of the initial Gaussian field, and if the power spectrum is smooth enough as CDM models. The 
genus asymmetry parameter for CDM models with shape parameter T, normalized by erg is listed for various smoothing 
length in Table 8. Since the factor - increase with the spectral index, as seen from Table 3, this factor 

also increase with the smoothing length for CDM model on scales of interest. On the other hand, the factor gq is a 
decreasing function of the smoothing length. As a result, the genus asymmetry parameter is not a monotonic function of 
the smoothing length. For example, the genus asymmetry parameter in the V — 0.2, erg = 1 CDM model is approximately 
0.08 on wide range of smoothing length of 10/i _1 Mpc ^ R 20/i _1 Mpc. 

From the second order formula for the genus of equation (3.53), the genus at the troughs is proportional to — 2±v / 3(S 1 ^ 2 ^ — 
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S^)ao- Therefore, the fraction of the deviation from the Gaussian prediction atv — ±\/3 is ±y/3(S^ — S^)a /2. When 
the genus asymmetry parameter is 0.08, this fraction is about ±7%. Thus, the perturbation theory predicts that the values 
of genus at positive and negative troughs differ 14% for the T = 0.3, a$ = 1 CDM model, and so on. This difference 
increases with the shape parameter V. The degree of deviation is qualitatively consistent with the numerical result in 
Fig. 9 of Canavezes et al. (1998), although quantitative comparison is still difficult because of the noise in the simulation. 

The statistics for the velocity field is plotted in Figure 2. Since the the sign of the skewness parameters is negative 
in this case, the weakly nonlinear evolution of statistics against the density threshold v indicates the sponge-like shift. 
In terms of the volume-fraction threshold v, on the other hand, meat-ball shifts are observed even for relatively bluer 
spectrum with n s < 0. 

The 2D projected galaxy statistics are dependent on the selection function of galaxies and cosmological models. As 
an example, we assume the APM luminosity function (Laveday et al. 1992) for galaxies with B band magnitude limit 
wi; m = 19. The differential number count dN/dz(z) for this sample is plotted in Figure 3. The resulting mean redshift 
is (z) — 0.12. The relation between the differential number count and the normalized mean number density n(x) in 
comoving coordinates are given by 

« [xM] = dz r ~,„., ■ (5-2) 



where 

H{z) = H V(l + z) 3 n + (1 + z) 2 (l - fi - Ao) + A , (5.3) 
X( z ) — [ Z T77~~J\j (5-4) 



H(z') 

are the Hubble parameter and the comoving distance at redshift z, respectively. Once the selection function n(x) is given 
by equation (5.2), the integration of equations (4.55)-(4.57), and the interpolation of tabulated values of in Table 6 
give the skewness parameters. For power-law power spectra, the skewness parameters are given by the simpler integration 
of equation (4.78) and the values of C (a) in Table 7. 

(2) (2) 

In Figure 4, 2-dimcnsional statistics, N±, G2, and V 2 arc plotted with assumed cosmological parameters, VLq = 0.3, 
Ao = 0.7, and APM luminosity function with limiting magnitude 19. The nonlinear parameter is assumed as <tq = 0.2. 
For the CDM model, the shape parameter is assumed as T = 0.25, and we take the smoothing angle as 9{ = 1°. With this 
smoothing angle the value <r = 0.2 corresponds to the normalization er 8 = 1.22. Basic features for these 2-dimensional 
statistics are the same as for 3-dimensional density field, except that the smoothing scale in CDM model we adopt 
corresponds to smaller scale than in the example of 3-dimcnsional density field of the Figure 1. 



6. CONCLUSIONS 

In this work, we comprehensively presented a basic formalism to treat the statistics of smoothed cosmic fields in per- 
turbation theory. This formalism provides a methodology on evaluating how various statistics deviate from the prediction 
of simple random Gaussian fields. As long as the non-Gaussianity is weak, the behavior of the statistics caused by non- 
Gaussianity is predicted by our formalism, which enable us to quantitatively compare the statistical quantities and the 
source of the non-Gaussianity. This method is considered as an extension of the Edgeworth expansion, which has been 
proven to be useful in various fields of research, as long as the non-Gaussianity is weak. In this paper, we derive useful 
formulas and relations focusing on application of the second-order perturbation theory to various cosmic fields. 

Several examples of statistics of cosmic fields in second-order perturbation theory are investigated in datail, including 
level-crossing statistics, 2D and 3D genus statistics, 2D extrema statistics, and the Minkowski functionals, which are 
extensively used in cosmology. More complicated statistics, such as 2D and 3D density peaks, can also be calculated, 
although they are more tedius. 



Table 8 

Genus asymmetry parameter (S (2) — 5 (0) )o- for a CDM model with shape parameter r. The power spectrum is 

NORMALIZED BY (78 = 1- THE PARAMETER E IS SET AS E — 3/7. 

Smoothing length R [/i _1 Mpc] 
T 6.0 8.0 10.0 12.0 14.0 16.0 18.0 20.0 



0.1 .. 


. -0.101 


-0.045 


-0.011 


0.011 


0.025 


0.034 


0.041 


0.045 


0.2 .. 


. 0.017 


0.054 


0.071 


0.078 


0.081 
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0.079 


0.076 


0.3 .. 


. 0.088 


0.107 


0.111 


0.108 


0.102 


0.096 


0.090 


0.084 


0.4 .. 


. 0.138 
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0.134 


0.123 


0.113 


0.102 


0.093 


0.085 


0.5 .. 
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0.165 


0.149 


0.132 


0.118 


0.105 


0.094 


0.085 
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Fig. 3. — The differential number count for the APM luminosity function with B-band limiting magnitude 19. The normalization is 
arbitrary. 
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A particular interest in cosmology of our method is in the application to the cosmic fields. Even if the cosmic field was 
random Gaussian at the initial stage, the gravitational evolution induces the non-Gaussianity. The gravitational instability 
is a well-defined process, so that we can evaluate the non-Gaussianity without any ambiguity when the evolution remains 
in the quasi-linear regime provided that the biasing from the non-gravitational process is simple enough on large scales. 
Therefore, we performed the pcrturbative analysis to obtain the necessary skewness parameters based on the gravitational 
instability theory. We considered the 3D density field, the 3D velocity field, the 2D projected density field. In the 
application of second order theory to various statistics of smoothed cosmic fields, three types of skewness parameters are 
commonly useful, i.e., S^°\ an S^ 2 \ Extensive calculations of these parameters for various cosmic fields are one of 
the new results of this paper. It would be true that other skewness parameters are needed when other complex statistics 
are considered. Such other parameters, if needed, are similarly calculated by the method we outlined in this paper. 

We find the lowest order deviations from the Gaussian predictions of various statistics of smoothed cosmic fields depend 
only on the differences of the skewness parameters when we use a threshold v, which is rescaled by a volume-fraction 
of the smoothed field. This rescaling makes the lowest deviations much smaller than in the case of direct threshold v. 
This is because three skewness parameters take similar values if it is arisen from the gravitational evolution. For the 
phcnomcnological hierarchical model, these parameters arc identical. In this case, the weakly nonlinear correction of the 
statistical quantities in terms of the volume-fraction threshold vanishes. When evaluated by the second-order perturbation 
theory of density fluctuations, those three types of skewness parameter are still close to each other. This fact explains 
the smallness of the deviations from Gaussian predictions of statistical quantities like genus, level-crossing, or Minkowski 
functionals when the rescaled threshold by volume-fraction is used. 

We discussed small, but detectable deviations from Gaussian prediction of the 3D genus curve against rescaled threshold 
in detail. In the framework of the second-order perturbation theory, a prominent deviation of the genus curve occurs at 
the two troughs of the curve. Relative depths of left and right troughs in the genus curve show slightly meat-ball shift 
for CDM-like models. We found the degree of this asymmetry is proportional to the combination (S^ — S^)ao- We 
call this factor as a genus asymmetry parameter, which we propose as theory-motivated parameter that characterize the 
asymmetry of the genus curve. The genus asymmetry parameter can observationally be obtained by fitting the genus 
curve by Hermite polynomials as equation (5.1). Qualitative comparison with the numerical simulations in literatures 
suggests that such asymmetry can actually be observable. We have not estimated over what dynamic range in nonlinearity 
parameter like ao, wo, or v<jq are each perturbative expressions valid. This estimation requires a systematic comparison 
with large N-body simulations, which is beyond the scope of this paper, and will be given in a subsequent paper of the 
series. 

In principle, any order in the perturbation theory can be calculated as further as one would like. Although the 
computation of the higher-order theory becomes more and more tedious, the necessity of the comparison with large-scale 
cosmological observations is a good reason to perform such computation as further as we can. One of the spectacular 
example of the detailed comparison between perturbation theory and observations is the fine structure constant in quantum 
electrodynamics (e.g., Kinoshita 1996). Our analysis in this paper will be extended to the third-order perturbation theory 
in a subsequent paper of the series. The present time is in an unique decade when the observations of cosmic fields are in 
unforeseen progress, like large-scale redshift surveys, detailed mapping of CMB fluctuations, gravitational lensing surveys, 




STATISTICS OF SMOOTHED COSMIC FIELDS 



27 



and so forth. Statistics of smoothed density field with higher-order perturbation theory will provide an unique method 
to analyze those high-precision data. The precision cosmology is undoubtedly providing clues to unlock the door to the 
origin of the universe. 

I would like to thank M. Kerscher and B. Jain for discussions. I wish to acknowledge support from JSPS Postdoctoral 
Fellowships for Research Abroad, and from the Ministry of Education, Culture, Sports, Science, and Technology, Grant- 
in- Aid for Encouragement of Young Scientists, 13740150, 2001. 

APPENDIX 
A. USEFUL GAUSSIAN INTEGRALS 

In this appendix, we give Gaussian integrals which are useful in this paper. In the following, H n is the Hermite 
polynomials, 

H n {u)=e^(-^Xe-^, (Al) 



(A2) 



and we further employ the notation 

Several Hermite polynomials are 

H (v) = l, H 1 {v)=v, H 2 {u) = v 2 - 1, H 3 (u) = v z - 3i/, 

Hi{v) = v A - 6is 2 + 3, H 5 (u) = v h - 10z/ 3 + Ibv. (A3) 
The Hermite polynomial at zero is given by 

h m / °' (n : odd) 

Un ^>- \ (-l)»/ 2 (n - 1)!! (I : even) 

ee h n . (A4) 

We generalize the above definition of h n to the case n < by interpreting (n — 1)!! as the appropriate gamma function 
so that (—1)!! = 1, (—3)!! = —1, etc. For example, h-2 = 1, h = 1, h 2 = — 1, hi = 3, and so forth. In this appendix, 
we give useful Gaussian averages (• • -) G of equation (2.20) for normalized cosmic field a defined by equation (2.5) and its 
spatial derivatives m — a.i. In the following, r\ represents any one of the components m. First, concerning a, 



/ d k 5(a — v) \ e 



\-^r^^(a)J G = ^H k+n M, (A5) 



/d k 8{a-v) v (n \ _ e~f 1 

/2tt 



\ — d^~ Hn{a) ' =~^ H k+n-M- (A6) 



Second, concerning r], in the notation of equation (A4), 

'd l \ V \\ 12 ( 



drf / G V 71" \Vda 
/d l S{r])\ 



(A7) 



hi-2, (A8) 



\ 8 V l I q V27rVVda y 



-i-i 



h t . (A9) 



B. limber's equation for bispectrum 

The correlation functions on projected sky is expressible by the 3-dimcnsional correlation functions. The explicit 
relation for the two-point correlation function is given by Limber's equation (Limber 1954). The Fourier-space version of 
Limber's equation is given by Kaiser (Kaiser 1998) and is somewhat simpler. His argument was generalized to higher-order 
correlation functions and their Fourier transforms (Scoccimarro, Zaldarriaga, & Hui 1999; Buchalter, Kamionkowski, & 
Jaffe 2000). Since the higher-order Limber's equation in Fourier space was discussed in the context of weak lensing field 
in literatures, here we review the derivation in a way more useful in this paper. Let 2D projected field / be the projection 
of a time-dependent 3D field F(x; r) along a light-cone: 

f(0) = J d x S 2 (x)q(x)F(x, 8S k ( x );t - x), (Bl) 

where q(x) is some radial weighting function and x is t ne radial comoving distance, and r is the conformal time at the 
observer. The past light-cone of the observer is specified by the equation x = r — r. The comoving angular distance 
Sk (x) is defined by equation (4.46). 
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In the following, we explicitly derive the relation for 3-point correlation function, and bispectrum. From the Limber's 
equation, power spectrum is already derived by Kaiser (1998): 

PfH = J dxSUx)q 2 (x)P F (^- y T - x), (B2) 

where Pf and Pp is the power spectrum of fields / and F, respectively. We generalize this equation to the one for 3-point 
statistics below. The generalization of the following derivation to higher-order statistics is straightforward. The angular 
3-point correlation function w^p of / is 

wf ) (0 1 ,e 2 ,d 3 ) = J d X iS 2 K {xi)q{xi) J d X2 S 2 K { X 2)q{x2) J d X3 S£(x3)q(X3) 

x (F (xuOiSk (xi) ; r - X) F ( X 2,e 2 S K ( X2 ) ; r - x) F ( X 3, 3 S K (xs) ; r - x)) 
dxSK 6 (x)q 3 (x) / dxidxiCF (xi,OiSk(xi);X2,02Sk(X2);X,Q3Sk(x); t o - x) , (B3) 



where (f(xi; ■ ■ ■ ; x 3 ; r) is the spatial 3-point correlation function of the field F with the 3-point configuration (x\; x 2 ; x 3 ) 
at conformal time r. According to the spirit of the Limber's equation, we assume that S^(x)q(x) 1S slowly varying 
compared to the scale of the fluctuations of interest and also that these fluctuations occur on a scale much smaller than 
the curvature scale. The equation (B3) is the generalization of the Limber's equation to higher-order correlation functions. 

Now we transform equation (3.18) to obtain the 2D bispectrum. We use the following convention of the Fourier 
transforms 

/M = J d 2 9f(6)e-*" e , (B4) 

F(k; T ) = j d 3 xF(x;T)e-* kx , (B5) 
the bispectrum Bf of the 2D field /, and Bp of the 3D field F are defined by 

(/(wi)/(w 2 )/(a>3)) = (2^) 2 <5 2 (u;i +^2 +a>3)B/(wi,W2,W3), (B6) 

/ F(k 1 ;T)F(k 2 ; r)F(k 3 ; t)\ = (2^) 3 ,5 3 (fe 1 + k 2 + k 3 )B F (k u k 2 , fc 3 ; t), (B7) 



where = \oJi\ and fcj = |fej|. The Dirac's delta function comes from the translational invariance of statistics. From the 
relations, 

B f (w u w 2 ,w 3 ) = J d 2 M 2 02«/(0i,02,O)e-^- 01 -^- 02 , (w 3 = |wi+u;2|), (B8) 

Cf (*i,*2,*3;r) = J ^^B F (k u k 2 ,k 3 -,T)e ik ^-* s)+^-(— 3), ( B9 ) 

the Fourier transform of the equation (B3) reduces to a simple equation, 

= J d X S 2 M q>(x)Bp(^ y ^;r - x). (BIO) 

This is a bispectrum version of the Limber's equation and the generalization of the Kaiser's equation for power spectrum 
(B2). 

C. SYMBOL INDEX 

In Table C9, the quantites used in this paper are listed. 
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Table C9 
Symbol Index 



Symbol 


Definition 


Equation 


\ 1 


Ensemble average 


(2 18) 




Ensemble average by Gaussian field 


(2.20) 


'-'3 


2D spectral parameter 


f4 55) 


n 


Time-dependent density parameter 


(4.53) 




Density parameter 


(4.53) 


a 


Normalized field of / 


(2.5) 


cx^j , etc. 


d 2 a/dxidxj = didja, etc. 


(2.6) 


7 


spectral parameter 


(2.38) 


5 


Mass density contrast 


(4.81) 


^3D 


3D density contrast 


(4.48 


6 


Galaxy number density contrast 


(4.81) 




= a.ij , Tensor notation 




7- ■ 


1 ilYlPflT POTTl l»l Y~i a t 1 On OT Cv 'A Tl f\ t ■ ■ 
J_JlIl<^dl LUIllUlIlOj L1U11 UI iJi cLLLvl 






— cv ■ A/pptor notation 

L>c ^ ^ V v^V~> 1 11 vy U CJUKjIVJ LI. 







Ancmlar coordinator in a small "natch of the- slcv 

JVllii LtlGll UlllCl LVjO 111 Clj OlllClll LyClLiVjll Ul L1H_j urv V 


(A 48) 


K 


Local convergence field 


(A 79) 


\ 


Timp-rlp'npnrlpnt scaled rosmolopnoal constant, 

1 1111*^ l_7 v^lll_l\_vll U O CI 1 \^ \_L l_. VJOlllUlVJ ZLl^CAtL L<VJ110 U CX11 U 


(A 53) 


An 

A 


paler! cosmolocnpal constant 

UvQilbQ LUD111U1U cL LKjGiL LUIIl) u till u 


(A 53) 




n rfit: n r~\l H rnr iron ctwC'C* 
J_ 111 L.t>llUl(_l UV Vcxl IctH^L- 






Tnrpsholrl hv vaIhttip fraction 

J_ 111 COllUlU UV V \_)± LilllC 11 OAj LHJ11 


C\ 4Q) 


n 

y 


3D pomovinp - rlonsitv field 


(A 47) 


Hez 


2D Weighted extrema 


(3.29) 


Hp 


2D projected density field 


(A 47) 


CO 


Variance of / 


(2.4) 




Spectral parameters 


(2.30), (2.31) 


TO 


Conformal time at the present 


(4.47) 


X 


Comoving distance 


(4.46) 




Volume of the unit ball in k dimensions 


(3.40) 




Vector of spatial derivatives of a 


(2.6) 


B(k 1 ,k 2 ,k 3 ) 


Bispectrum at the present time 


(4.5) 


-B 2 d(wi,w 2 , w 3 ) 


2D bispectrum 


(4.50) 


-B3D(fci,fc2,fc 3 ;r) 


3D bispectrum at conformal time r 


(4.50) 


C(«) 


Integrand of the 2D skewness parameters 


(4.58) 




Integrand of the 2D skewness parameters 


(4.59) 




Contributions to the 2D skewness parameters 


(A 69) 




linear prowth rate 

11A1^_-C11 tl WW Lill IClUC 


(A 51) 




Parameter of the nortnrhation thoorv 

_1 Cll ClrlllV_y LO-L VJ1 ullL VJ \j 1 Li HI UCX L1U11 L 11L.L/1 y 


(A 12) 


E v 


Parameter of the perturbation theory 


(4.38) 


F 


Cosmic field, in general 


(2.18) 




d^F/dA^dA^dA^, etc. 


(2.22) 


F(a,/3,7;z) 


Gauss hypergeometric function 2 f\ 


(4.30) 


G 2 


2D genus statistic 


(3.13), (3.14) 


G 3 


3D genus statistic 


(3.20), (3.22) 


i? 


Time-dependent Hubble parameter 


(4.36) 


H n 
"0 


nuuuie consiaiiL 




/„ 


Modified Bessel function 


(4.17) 




Fourier counterpart of A^ 


(2.7) 


Jl 2) ({m i3 }) 


Integers defined by Table 1 


(3.24) 


K 


Spatial curvature 


(4.46) 


H n {v) 


Hermite polynomials 


(Al) 


H-M 


Extended Hermite polynomial of order —1 


(A2) 


M {n) 


n-th order cumulants of A^ 


(2.8)-(2.12) 


M (n) 


Normalized cumulants of A^ 


(2.21) 
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Table C9 — Continued 



Symbol 


Definition 


Equation 


A/f 

1V± 


qooc\v\ r\ (~\yr\ ot on m i i 1 a nt~ r^T /a 

<J± LIL-l L-LlIllUlclIlL Ul j!^ 




iV i 


T ^M/TM-PTTlQQin O" QtfltlQtip 

i_ji_.vli i_.i usaing & Lxabio 




l v 2 


ULllg Ull OLdrLlBLlL 


CS 4) 


1 v 3 


Area statistic 




Pn 


A/n i 1 tiva n a tp (TfinQt;i9Ti ni Qtri ni i t i nn function 

±V1 111 LI V OjI IdlL VJTOj LloolClll UIO LI 1UUL1U11 1 l_J.li.l_- Ll^Jll 


(2 7) 


P 


A/fnltivnriatc "nrot>?) t>litiv distribution function 

1V1 U.1 Lil V Clrl 1 Cli L yjL\J l_/Clrl_/ll LjI y UlO Lilll_/HLilWll 1 Llll"^ UJU11 


(2 7) 




Power Qnpptrnm at rtrcQent 

J. <_JW<^1 OlJ^L-Ll Hill CLL JJ1 L^iZl^ll L 


(4 4) 


-» 2D ^ ,1 




(4 4Q\ 




Xl) nrra^PT Qncctrnm at frminrmfll time t 

»J 1 7 l.y*_JVV'^l kJ I J Kj L 1 Hill CI b L/U111U1 llltll Lj111H_/ / 


(4 4()) 


L1JN V A v 


Tiincpir nnwpr miprtrnm 

J_ilAl^CAjl UL/VV^l GL/L.ljlil Utlll 


(4 10i 




TTicrarcnical arrrnlitnrlc of trie /{-"point function 


(4 li 




Skewness parameters 


(2.61)-(2.64) 




Integrand of the skewness parameters 


(4.15) 


qa/3 ( T> \ 


Contributions to the skewness parameters 


(A 91 'I 




Galaxy skewness parameters 


(4.46) 


S K 


Comoving angular diameter distance 


(4.46) 


TcDM 


CDM transfer function 


(4.25) 




Minkowski functionals 


(3.32), (3.33), (3.37), (3.40) 


w 


Smoothing function in real space 


(2.1) 


W R 


Smoothing function in real space 


(2.1) 


z 


Generating function 


(2.7) 


b 


Bias parameter 


(4.81) 




Nonlinear bias parameter 


(4.81) 


d 


Dimension of the sample space 


(2.1) 


dN/dz 


Differential number count 


(5.2) 


I 


Smoothed cosmic field, in general 


(2.1) 


fv 


Volume fraction on the high-density side 


(3.49) 


9f 


Logarithmic derivative of the growth factor 


(4.39) 


hi 


Hermite polynomial at zero, Hi(0), extended to I < 


(A4) 


n 


Selection function per unit comoving volume 


(4.47) 


n s 


index for the power-law power spectrum 


(4.26) 


x, y, z 


Linear combinations of a and Qj's 


(2.39)-(2.41) 
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